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Sherwood.  Theoretical  Meeting 
March  22-23,  1971 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 


All  sessions  will  be  held  in  Room  109  Warren  Weaver  Hall 
[251  Mercer  Street,  4th  Street  East  of  Washington  Square] 
except  session  IV-b  in  Room  102. 

Where  known,  underlined  author  will  present  paper. 
Opening,  9  ajn>  Monday,  March  22 

Welcome;  Louis  Nirenberg,  Director,  Courant  Institute 

Announcements 
Session  I,  Chairman  -  C.  L.  Hedrick 

1.  W.  A.  Newcomb.  "Stationary  Equilibria  of  a  Guiding- 
Center  Plasma'  (15  min. ) 

2.  D.  Lortz,  E.  Rebhan,  G.  Spies,  "Sufficient  Criteria 
for  Magnetohydrodynamic  Stability"  (15  min.) 

3a.  R.  Ga.iewski,  "mhD  Equilibrium  of  an  Elliptical  Plasma 
Cylinder 

b.  R.  Dagazian,  B.  Coppi,  "MHD  Instabilities  in  Confinement 
Configurations  with  Shaped  Magnetic  Surfaces"  (20  min.) 

4.  N.  Friedman,  "Free  Boundary  Scylla  Equilibria"  (15  min.) 
COFFEE 

5.  H.  Grad,  A.  Kadish,  D.  Stevens,  "Free  Boundary  Tokomak 
Equilibria  with  No  Outer  Wall"  (15  min.) 

6.  J.  Freidberg,  "Stability  of  the  Straight  £  =  1   Scyllac 
Configuration"  (15  min.) 

7.  H.  Grad,  J.  Marsh,  H.  Weitzner,  "A  Spectral  Analysis 
of  the  Screw  Pinch"  (15  min. ) 

Rl.   Report  from  Princeton,  J.  Dawson,  P.  Rutherford  (30  min.) 

Session  II,  2  pm.  Chairman  -  J.  Freidberg 

R2.   Report  from  Los  Alamos,  R.  Morse  (20  min.) 

8.  J.  D.  Callen,  "Tokomak  Kinetic  Theory;  Trapped- Particle 
Modes"  (15  min. ) 
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Session  II  (cont'd) 

9.   A.,  k.    War.e,  A.  B.  Macmahon,  "Trapped  Particle  Contri- 
butions to  the  Heating  in  a  Tokomak"  (15  min. ) 

10.  3.  Yoshikawa,  J.  Schmidt,  "Anomalous  Viscosity  as  a 
Possible  Explanation  for  an  Anomalous  Skin  Effect   (15  min.) 

COFFEE 

11.  F.  L.  Hlnton,  C.  W.  Horton,  Jr.,  "Bootstrap  Current- 
Driven  Drift  Instability  in  Tokomaks"  (15  min.) 


12. 


D.  K.  Bhadra,  C.  S.  Liu,  T.  Ohkawa,  "Neoclassical 
Diffusion  in  Axisymraetric  Torus  in  the  "intermediate 
Region""  (15  min.) 


13.  N.  K.  Winsor,  E.  C.  Bowers,  M.  A.  Hellberg,  J.  M.  Dawson, 
"Shocks  and  Diffusion  in  Low-P  Toroidal  Plasmas"  (15  min.) 

14.  W.  Ross,  M.  N.  Rosenbluth,  "Stabilization  of  Dissipative 
Trapped  Particle  Instability"  (15  min.) 

R3.   Report  from  N.  Y.  U. ,  H.  Grad  (20  min.) 

Session  III,  9  am,  Tuesday  March  23,  Chairman  -  N.  Winsor 

r4.   Report  from  M.  I.  T.,  B.  Coppi  (20  min.) 

15.  D.  V.  Anderson,  J.  Killeen,  "Computation  of  Finite-Beta 
Plasma  Equilibria  in  Minimum-B  Mirror  Systems"  (15  min.) 

16.  D.  Dobrott,  H.  Grad,  "Llmiter  Induced  Perturbations 
and  the  Small  Gyro  Radius  Limit"  (15  min.) 

17a.  B.  D.  Fried.  C.  F.  Kennel,  "Low  Frequency  Drift  Cyclotron 
Instability  ' 

b.  D.  W.  Forslund,  R.  L.  Morse.  C.  W.  Nielson,  "Electron 
Cyclotron  Drift  Instability^'  (20  min.) 

COFFEE 

18.  E.  J.  Valeo,  C.  Oberman,  W.  L.  Kruer,  "Nonlinear  Theory 
for  Plasma  Driven  by  a  Large-Amplitude,  High-Frequency 
Field"  (15  min.) 

19.  C.  L.  Hedrlck,  Jr.  "Suppression  of  the  I\niistler  Instability 
by  Relativistic  Effects   (15  min.) 

20.  J.  P.  Freidberg,  B.  M.  Marder,  "High  Frequency  A.  C. 
Electrostatic  Plasma  Instabilities   (15  min.) 

21.  P.  K.  Kaw,  Y.  C.  Lee,  "Dynamic  Stabilization  of  Drift 

Waves  by  Radial  and  Azimuthal  High-Frequency  Fields"  (15  min.) 

R5.   Report  from  Oak  Ridge,  G.  Guest  (20  min.) 


■IV- 


Session  IV-a,  2  pm.  Chairman  -  D.  Dobrott 

R6.   Report  from  LRL,  D.  Pearlstein  (20  mln. ) 
R7-   Report  from  Texas,  W.  E.  Drummond  (20  min. ) 

22.  H.  R.  Strauss,  "High  Pressure  Tokomak  Equilibriuun"  (15  min.) 

23.  J.  M.  Dawson,  H.  P.  Furth,  F.  H.  Tenney,  "Thermonuclear 
Power  Production  by  Nonmaxwellian  Ions  in  a  Closed 
Magnetic  Field  Configuration"  (15  min.) 

24.  M.  S.  Chu,  "Thermonuclear  Reaction  Waves  at  High 
Densities'  (15  min.) 

COFFEE  3:^5  pm 

25.  J.  P.  Holdren,  A.  H.  Futch,  Jr.,  J.  Killeen,  A.  A.  Mirin 
R.  F.  Post,  "Multi-Species  Fokker  Planck  Calculations 
and  Q  Values  for  D-T  and  D-He  Mirror  Reactors"  (15  min.) 


26.   J.  M.  Kindel,  F.  W.  Perkins,  "Will  Ion  Cyclotron  Waves 
Heat  Tokomak  Plasmas?"  (15  min.) 

27a.  R.  Lee,  R.  N.  Sudan,  "Cross  Field  Injection  of  a 

Relativistic  Electron  Beam  into  a  Magnetized  Plasma" 

b.  R.  V.  Lovelace,.  R.  N.  Sudan,  "Precession  of  a  Strong 
E- Layer"  (20  minT) 

23.   N.  C.  Chris tofilos,  "E-Layer  Equilibria  in  the  Presence 
of  a  Toroidal  Field"  (15  min.) 

Session  IV-b,  2:45  pm.  Room  102,  Chairman  -  R.  Gajewski 

29.  P.  H.  Sakanaka.  T.  C.  Marshall,  C.  K.  Chu,  "Formation  and 
Interaction  of  Ion  Acoustic  Solitary  Waves  in  a 
Collisionless  Warm  Plasma"  (12  min.) 

30.  P.  N.  Hu,  "Evolution  of  Shock  and  Nonlinear  Compression 
Waves"  (12  rain. ) 

31.  E.  A.  Williams,  C.  Oberman,  "Noise  Broadening  of  Trapped- 
Particle  Echoes"  (12  min.) 

32.  T.  M.  O'Neil,  J.  H.  Winfrey,  "Nonlinear  Interaction  of  a 
Small  Cold  Beam  and  a  Plasma,  Part  II  (Random  Phase 
Approximation  Revisited)"  (12  min.) 

COFFEE  3:45  pm 

33.  M.  Bineau,  "On  The  Existence  of  Force  Free  Magnetic 
Fields  in  a  Simply  Connected  Domain"  (12  min.) 

34.  R.  E.  Aamodt,  "Containment  in  Open-Ended  Systems  in  the 
Presence  of  Fluctuating  Electric  Fields"  (12  rain.) 
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Session  IV-b  (cont'd) 

35.  M.  N.  Rosenbluth,  P,  H.  Rutherford,  R.„Hazeltlne, 
"Effect  of  Field  Asjoranetries  on  Neoclassical  Confinement 
in  Tokomaks"  (12  min. ) 

36.  R.  L.  Dewar  "A  Lagranglan  Theory  for  Nonlinear  Wave- 
packets  in  a  Colllsionless  Plasma"  (12  min.) 

37.  F.  V.  Coroniti,  B.  D.  Fried,  R.  B.  White,  "Shock  Wave 
Structure  in  a  Two  Component  Plasma"  (12  min.) 

Supplementary  Program  (by  Title  or  in  Review  Papers) 

C.  D.  Striffler,  T.  Kammash,  "Evolution  of  Unstable 
Bernstein  Modes  in  Astron  Plasma" 

R.  K.  Varma,  C.  W.  Horton,  Jr.,  "A  Schrodinger  Description  of 
Non-Adiabatic  Particle  Loss  and  of  Magnetic  Surfaces" 

J.  K.  Percus,  G.  J.  Yevick,  "Decay  of  Pair  Correlations  in  a 
Classical  Fluid" 

U.  Daybelge,  "Neoclassical  Electrical  Conductivity  of  Plasmas 
with  Trapped  Particles" 

G.  Schmidt,  "Self  Consistent  Field  Theory  of  Relativistic 
Electron  Rings" 

F.  Winterberg,  "Magnetically  Insulated  Ultra-Hi^h  Voltage 

Transformer  and  Controlled  Fusion  Research 

S.  H.  Schneider,  C.  K.  Chu,  "Numerical  Simulation  of  Strong 
Plasma  Shock  Waves  Produced  in  an  Electromagnetic  Shock 
Tube" 

B.  Rosen,  G.  Schmidt,  "The  Role  of  Trapped  Particles  in 
Plasma  Instabilities" 

W.  B.  Thompson,  "Low  Frequency  Heating" 

M.  Lampe,  J.  B.  McBride,  W.  M.  Manheimer,  J.  H.  Orens,  R.  N.  Sudan, 
"Theory  and  Simulation  of  the  Beam  Cyclotron  Instability" 

W.  L.  Kruer,  J.  M.  Dawson,  "Electron  Vortices  and  Anomalous 
DC  Resistivity" 

P.  K.  Kaw,  R.  M.  Kulsrud,  "Relativistic  Particle  Motion  in 
Super- Intense  Laser  Beams" 

E.  A.  Frieman,  J.  M.  Greene,  K.  E.  Weimer,  "Toroidal  Effects  on 
Kink  Modes  in  Tokomaks 

H.  Okuda,  B.  Rosen,  J.  M.  Dawson,  "Computer  Simulation  of 
Drift-Cyclotron  Instability^' 

J.  T.  Yen,  "Analysis  of  Nonlinear  Plasma  Waves" 

G.  Bateman,  "Energy  Principle  with  Given  Distributed  Inductance" 

T.  Yeh,  "Hydromagnetic  Stability  of  Bifurcated  Equilibria" 
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stationary  Equilibria  of  a  Guiding-Center  Plasma 
William  A.  Newcomb 
University  of  California,  Lawrence  Radiation  Laboratory 

Livermore,  California 

We  consider  equilibrium  states  of  a  guiding-center  plasma, 
assuming  either  _v  =  0  (special  case)  or  _v;  ?^  0  (general  case), 
where  v;  is  the  drift  velocity  (=  E  x  B/B  )  .   Let  H  be  the  total 
energy.   Assuming  _v  =  0,  the  variational  condition  for 
equilibrium  may  be  given  as  5H  =  0,  provided  that  5J  =  0  is 
imposed  as  a  side  condition,  where  J  is  the  areal  invariant  in 
the  two-dimensional  phase  plane  of  the  particle  motion  in  the 
parallel  direction.   This  was  first  shown  by  Grad,  and  we 
generalize  it  now  for  nonvanishing  v.   In  the  general  case, 
the  equilibrium  condition  is  again  5H  =  0,  but  now  with  two 
side  conditions:   6J  =  0  (as  before),  and  6C  =  0,  where  C  is  a 
modified  circulation  integral  for  the  macroscopic  _v  motion. 
The  explicit  form  of  C  is  derived  from  a  symmetry  principle: 
invariance  under  excahnge  of  equivalent  flux  lines.   (In  the 
same  way,  the  areal  invariant  arises  from  the  principle  of 
invariance  under  exchange  of  equivalent  particles,  so  that  a 
unified  treatment  of  both  the  side  conditions  may  be  given.) 

-X- 

Work  performed  under  the  auspices  of  the  U.S.  Atomic  Energy 
Commission. 
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Sufficiclent  Criteria  for  Magnetohydrodynamic  Stability 
D.  Lortz,  E.  Rebhan,  G.  Spies 
Max- Planck- Institut  fur  Plasmaphysik,  Euratom  Association 
Garching,  Fed.  Rep.  of  Germany 

A  class  of  sufficient  magnetohydrodynamic  stability  criteria 
has  been  derived  for  equilibria  which  are  periodic  along  the 
magnetic  field  lines  (i.e.  having  some  symmetry  or  closed  field 
lines) .   The  sufficient  criteria  known  hitherto  are  contained 
as  special  cases.   In  this  class  there  is  an  optimum  criterion 
which,  in  general,  is  less  restrictive  than  the  former  sufficient 
criteria.   However,  its  application  requires  solving  an  eigenvalue 
problem.   More  convenient  criteria  have  been  derived  from  this 
optimum  criterion  for  helically  symmetric  systems,  including  the 
limiting  cases  of  axial  and  plane  symmetry.   These  criteria 
appear  well  suited  to  calculating  lower  bounds  for  a  critical  P. 
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MHD  Equilibrium  of  an  Elliptical  Plasma  Cylinder 
Ryszard  Gajewski 
Departmeat  of  Physics  and  Research  Laboratory  of  Electronics,  M.I.T. 

Cambridge,  Massachusetts  02139 

It  has  been  suggested  that  departure  from  the  circular  shape 
of  the  toroidal  plasma  column  cross  section  may  be  advantageous  to 
stability.   One  way  of  achieving  the  desired  shape  of  the  column 
cross  section  is  to  properly  shape  the  conducting  shell  which 
must  coincide  with  one  of  the  magnetic  surfaces  in  vacuum. 

In  the  present  study  magnetic  surfaces  are  found  analytically 
in  the  vacuum  region  surrounding  a  cylindrical  plasma  column  at  MHD 
equilibrium.   The  cross  section  of  the  column  is  assumed  an  ellipse 
and  the  longitudinal  current  density  is  assumed  constant  over  the 
cross  section.   The  magnetic  surfaces  are  deformed  ellipses,  elongated 
along  the  major  axis  (Fig.  1),  and  bound  by  a  separatrix  (Fig.  2). 
The  separatrix  tends  to  infinity  when  the  ellipticity  e  -►  0,  and 
approaches  the  surface  of  the  column  as  e  -♦  1.   The  case  when  the 
longitudinal  current  is  confined  to  the  surface  of  the  coliamn  is  also 
discussed. 
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Fig.  1 


Fig.  2 


This  work  was  supported  by  the  U.  3.  Atomic  Energy  Commission 
(Contract  AT(30-l)-3980) . 

'"L.  3.  Solovev,  V.  D.  Shafranov,  E.  I.  Yurchenko,  Nucl.  Fusion, 
Special  Supplement  1969,  P-  '55. 
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MHD  Instabilities  in  Confinement  Configurations 
with  Shaped  Magnetic  Surfaces 

R.  Dagazian  and  B.  Coppi 

Massachusetts  Institute  of  Technology 

The  equation  of  motion  is  solved  for  a  cylinder  of 
plasma  with  an  elliptically  or  triangularly  deformed  cross 
section  immersed  In  a  strong  magnetic  field.   A  dispersion 
relation  is  obtained  for  "kink"  like  perturbations.   In 
particular,  the  stability  of  the  plasma  column  under  such 
perturbations  is  studied  in  the  presence  of  the  metallic 
casing  which  serves  to  shape  the  magnetic  surfaces. 

* 
Work  supported  in  part  by  the  U.S.  Atomic  Energy  Commission, 

Contract  AT(30-1 )-598o. 


-4- 


Free  Boundary  Scylla  Equilibria 
N,  Friedman 
Cornell  University 
Ithaca,  New  York 

Straight  helical  equilibria  have  been  computed  using  a 
free  boundary  flat  pressure  model.   Given  the  shape  of  the 
outer  fixed  conductor,  the  shape  of  the  plasma  is  calculated. 
One  special  case  is  i  =  1,  which  is  pertinent  to  the  Scylla 
device.   In  this  case,  the  outer  boundary  is  specified  to  be 
an  eccentric  circle,  and  this  eccentricity  determines  the 
helical  field  perturbation.   Additional  parameters  are  the 
helical  wavelength,  the  ratio  of  the  radius  of  the  plasma  to 
that  of  the  wall,  and  the  magnitude  of  the  longitudinal  (pinch) 
current.   The  plasma  shape  and  displacement  are  calculated  as  a 
function  of  these  parameters.   The  eccentricity  creates  a  mutual 
Inductance  and  from  this  we  can  compute  either  the  strength  of 
the  induced  z-current  if  the  azimuthal  flux  is  kept  fixed  or  the 
induced  axial  EMF  if  the  z-current  is  held  at  zero. 
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Free  Boundary  Tokomak  Equilibria  with  No  Outer  Wall 
Harold  Grad,  Abraham  Kadish,  Donald  Stevens 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

A  two  dimensional  idealization  for  perfectly  conducting 
Tokomak  equilibria  has  been  studied.   Containment  is  provided 
by  the  vertical  field  alone;  there  are  no  outer  conducting 
walls.   Except  for  trivial  scaling,  there  is  a  one  parameter 
family  of  free  boundary  equilibria  specified  by  the  aspect 
ratio  or  by  the  ratio  of  vertical  field  strength  to  the 
constant  value  of  the  field  at  the  plasma.   For  small  vertical 
field,  the  aspect  ratio  is  large,  and  the  plasma  cross-section 
is  approximately  circular.   With  Increasing  vertical  field, 
the  plasma  shape  distorts,  and  the  axial  plasma  dimension 
becomes  much  larger  than  the  radial  dimension.   The  results 
are  obtained  by  conformal  mapping  and  are  explicit  in  terms  of 
elliptic  integrals.   Asymptotic  formulas  have  been  obtained 
for  large  and  small  aspect  ratio  and  numerical  values  in  between. 
Curves  are  given  for  the  variation  of  all  plasma  dimensions, 
temperature,  pressure,  density,  and  current  as  the  vertical 
field  is  varied  adlabatically  (as  in  the  proposed  ATC  experiment), 

These  results  are  preliminary  to  the  calculation  of  actual 
toroidal  free  boundary  equilibria  without  walls. 


stability  of  the  Straight  &   =  1   Scyllac  Configuration 

J.  P.  Freidberg 
Los  Alamos  Scientific  Laboratory 

The  stability  of  a  straight  6   pinch  with  a  superimposed  i   =  1 
field  is  studied  for  the  case  of  arbitrary  P,  in  the  sharp  boundary 
model.   Previous  investigations  of  this  problem  have  relied  on  one 
of  two  expansions  depending  upon  the  relative  size  of  two  dimen- 
sionless  parameters,  e  =  ha  s  pitch  number  x  plasma  radius  and 
6  =  shift/plasma  radius.   In  the  first  expansion  it  is  assumed  that 
S  <<  e  <  1.   For  this  case,  leading  order  terms  in  o  and  leading 
and  first  order  terms  in  e  are  kept.   The  plasma  is  found  to  be 
weakly  unstable  to  the  m  =  1  mode.   However,  substituting  into 
the  toroidal  equilibrium  relations  for  Scyllac  parameters  results 
in  a  violation  of  the  original  ordering;  that  is  for  equilibrium 
it  is  necessary  that  5-1,  e  ~  .1.   This  led  to  the  development 
of  the  second  expansion  in  which  it  is  assumed  that  e  <<  S  <  1. 
In  this  expansion  leading  order  terms  in  £  and  leading  and  first 
order  terms  in  5  are  kept.   The  results  of  this  theory  also  predicted 
weak  instability  although  for  typical  parameters,  the  growth  rate  was 
considerably  smaller  than  that  of  the  first  expansion.   The  calcu- 
lation presented  here  is  based  on  the  second  expansion  which  is 
more  closely  related  to  experimental  conditions.   By  calculating 
growth  rates  numerically  we  are  able  to  carry  the  expansion  to 
leading  and  first  order  terms  in  e  and  arbitrary  6.   The  main  result 
of  the  calculation  is  that  for  experimental  parameters  the  additional 
corrections  to  the  second  expansion  appear  to  be  very  important  and 
somewhat  suprisingly  predict  growth  rates  very  close  to  those  of  the 
first  expansion  except  for  very  high  P. 

Work  performed  under  the  auspices  of  the  U.  S.  Atomic  Energy 
Com-nission. 
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A  Spectral  Analysis  of  the  Screw  Pinch 
Harold  Grad,  Jeffrey  Marsh  and  Harold  Weitzner 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

Previous  analyses  of  the  stability  of  the  screw  pinch 
equilibrium  in  ideal  MHD  theory  have  been  done  in  many 
different  ways  and  are  the  source  of  many  qualitative  concepts. 
In  order  to  learn  more  about  this  simple  system  we  analyse  the 
stability  directly  from  the  equations  of  motion  and  determine 
the  spectrum  of  oscillation  and  growth  frequencies. 

Typical  of  far  more  general  topological  toroidal  con- 
figurations we  find  first  that  there  are  instabilities  not 
observable  in  the  usual  5w  formulation.   While  such  instabilities 
grow  algebraically  with  time  and  not  exponentially  they  may  well 
be  of  some  significance.   Even  the  ordinary  stable  theta  pinch 
may  exhibit  such  instabilities.   A  nonlinear  study  of  such 
instabilities  is  also  given. 

For  a  given  k  and  m  number  of  the  perturbation  the  continuous 
and  stable  and  unstable  point  spectra  are  described,  the  role  of 
the  Suydam  criterion  in  terms  of  the  spectrum  is  given.   For 
systems  only  slightly  Suydam  unstable  growth  rates  associated 
with  the  Suydam  modes  are  estimated.   Various  pathological  cases 
are  explored  for  significance.   The  applicibility  of  perturbation 
theory  to  various  instabilities  is  explored.   Possible  relationships 
between  mathematical  pathology  and  physical  significance  are 
considered. 
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Tokamak  Kinetic  Theory;  Trapped-Partlcle  Modes 

J.  D.  Callen 

Department  of  Aeronautics  and  Astronautics 

Massachusetts  Institute  of  Technology,  Cambridge,  Mass. 

An  Integral  equation  governing  kinetic  modes  in  a  colllslonless, 
current-carrying,  Tokamak-type  plasma,  which  Includes  magnetic  shear 
and  finite  "banana-width"  effects,  has  been  derived.   The  derivation 
is  accomplished  by  using  the  magnetic  flux  coordinates  of  an  axi- 
symmetric  system.   Namely,  we  use  f,  9,  ^  in  which  ip   is  the  poloidal 
flux  function,  C  the  toroidal  (axisymmetry)  angle  and  9   a  poloidal 
angle  defined  in  terms  of  the  surfaces  (x)  orthogonal  to  the  Tp 
surfaces  in  such  a  way  that  B  =  V^  x  V(q9-  C)  where  q  =  q(^)  =  2Tr/L 
is  the  "safety  factor".   The  particle  orbits,  including  the  gyro-, 
bounce  and  drift  motions,  are  derived  in  terms  of  these  coordinates. 
Integrating  the  perturbed  Vlasov  equation  along  these  orbits  for 
perturbations  of  the  type  *  .(^,  9)exp(i^C  -  inS  -  icot)  and  substituting 
the  result  into  Poisson's  equation,  we  obtain  the  integro-differential 
equation  for  0.(^,0)  governing  electrostatic  perturbations  in  an 

ajcisymmetric  system.   The  equation  is  similar  to  those  obtained  for 

1  2 

axisymmetric  multipoles   and  axisymraetric  magnetic  mirror  systems. 

It  differs  from  those  works  in  that:   1)  magnetic  shear  is  explicitly 

included;  2)  finite  gyroradius  and  finite  "banana-width"  corrections 

* 
Research  supported  by  the  National  Science  Foundation  (Grant 

GP-9557)  and  in  part  by  the  U.  S.  Atomic  Energy  Commission 

(AT(30-l)-3980). 

■"■p.  H.  Rutherford  and  E.  A.  Frieman,  Phys .  Flu.  11,  569  (1968); 
J.  B.  Taylor  and  R.  J.  Hastie,  Plasma  Phys.  10,  ^79  (1968). 

^C.  ¥.  Horton,  J.  D.  Callen  and  M.  N.  Rosenbluth  (to  be  published). 
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in  the  i'   motion  are  kept;  and  3)  the  magnetic  field  is  not  a 

vacu'jm  one,    but  rather  one  which  allows  for  a  plasma  current. 

We  have  used  this  new  integral  equation  to  investigate 

the  trapped-particle  interchange  modes-^  (<d  <<  ^^>    the  ion  bounce 

frequency).   In  particular,  we  study  the  radial  (or  Tp)    dependence 

of  the  unstable  modes  and  the  effects  of  shear,  removing  the 

ox  . 
previous  assumptions  that  -^  <<  |iq-n|  «  1  (close,  but  not  on 

a  mode-rational  surface).   The  f,9   eigenfunction  problem  is  solved 

iteratively  by  first  solving  the  9  (integral)  equation  and  then 

solving  the  ^  (differential)  equation.   Generally  speaking,  the 

unstable  modes  are  found  to  follow  the  magnetic  field  structure 

in  such  a  way  that  they  are  always  as  flute-like  as  possible. 

High  S,   number  modes  are  found  to  be  stabilized  by  the  finite 

banana-width  effects. 


^3.  B.  Kadomtsev  and  0.  P.  Pogutse,  Soviet  Phys.  -  JETP  24, 
1172  (1967). 
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Trapped  Particle  Contributions  to  the  Heating  In  a  Tokamak. 
A.  A,  Ware  and  A.  B.  Macmahon 
University  of  Texas  at  Austin 

Using  the  transport  equations  of  Rosenbluth  and  Hlnton, 
the  steady  state  solutions  have  been  obtained  for  n,  T  ,  T.  and 
B„  allowing  for  ohmlc  heating  and  the  various  trapped  particle 
effects.   Assuming  the  aspect  ratio  Is  fixed  and  utilizing  a 
dlmenslonless  form  for  the  equations,  the  single  parameter  which 
determines  the  character  of  the  solutions  is  a"  ./c,  where  a  is 
the  minor  radius  and  ^  .  the  ion  plasma  frequency;  this  parameter 
controls  the  degree  of  thermal  coupling  between  electrons  and 
ions.   The  properties  of  these  solutions  will  be  discussed. 

In  the  Tokamak  experiments  T-3  and  ST,  consideration  of 
the  time  constants  Involved  indicate  that  their  pulse  lengths 
are  too  short  for  such  a  steady  state  to  be  attained.   Compression 
and  heating  continue  throughout  the  pulse.   Initially,  ohmlc 
heating  will  be  the  dominant  heating  process,  but  once  P^  becomes 
comparable  with  unity  the  dominant  process  will  be  compressional 
heating  caused  by  the  trapped  particle  pinch  effect.   Approximate 
calculations  show  that  a  sharp  spatial  peak  in  temperature  can 
result,  as  observed  in  ST.   The  plasma  in  T-3  does  not  reach  this 
stage  because  of  the  longer  time  scale  and  the  fact  that  the 
driving  force  for  the  compression  (the  external  electric  field)  is 
reversed. 
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Anomalous  Viscosity  as  a  Possible  Explanation  for  an  Anomalous  Skin  Effect 

S.  Yoshikawa  and  J.  Schmidt 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  New  Jersey 

Anomalous  viscosity  of  electrons  across  the  magnetic  field 
may  be  responsible  for  the  observed  anomalous  skin  effect  in 
tokamaks .   Here,  a  model  is  proposed  to  determine  the  anomalous 
viscosity.   This  may  have  relevance  in  collisionless  shocks  across 
the  magnetic  field.   In  the  model,  electrons  which  drift  parallel 
to  B,  have  a  velocity  shear  in  the  perpendicular  direction  (say  x) . 
The  electrons  excite  waves  which  have  the  group  velocity  finite 
in  the  x  direction.   Then  the  waves  which  get  the  momentum  from 
the  electron  drift  propagate  to  the  region  where  the  drift  velocity 
of  electrons  is  smaller  than  the  phase  velocity  of  waves  parallel 
to  B.   Then  the  momentum  is  transferred  to  those  slow  electrons. 
That  is,  waves  act  as  carriers  to  transfer  the  momentum  between 
spatially-separated  electrons.   In  practice,  the  mean  free  path 
of  the 'wave  perpendicular  to  B,  i,  is  estimated.   The  momentum 
transfer  rate  between  electrons  and  waves,  r  ,  also  is  calculated 

from  quasilinear  theory.   The  anomalous  (kinematic)  viscosity  is 

2  -1 
then  estimated  as  ^  t   .   The  amplitude  of  the  wave  is  left  as  a 

parameter  to  be  determined  by  full  nonlinear  theory  and/or  by 

experiments.   This  work  was  supported  by  U.S.  Atomic  Energy 

Commission  Contract  AT( 30-l)-1238. 
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Bootstrap  Current-Driven  Drift  Instability  in  Tokamaks 
F.  L.  Hinton  and  C.  W.  Horton,  Jr. 
University  of  Texas  at  Austin 

The  stability  of  collisionless  drift  waves  is  examined  for 
Tokamak  geometry.   The  equilibrium  distribution  functions  are 
those  given  by  neoclassical  transport  theory  in  the  banana  regime. 
These  distribution  functions  are  known  to  give  a  "bootstrap" 
current,  carried  by  trapped  particles,  due  to  the  spatially  non- 
uniform distribution  of  bananas.   This  current  is  shown  to  be  the 
dominant  destabilizing  effect  for  drift  waves  (o)  jn;  cu^)  with  m- 
numbers  in  the  range  1  <<  m  <  B Vb„,  where  m  is  the  quantum 
number  of  the  mode  variation  that  short  way  around  the  torus. 
If  the  electron  to  ion  temperature  ratio  satisfies  yi  /T.'  >  r  /a.g, 
where  r  is  the  density  scale  length  and  a. g  is  the  average  ion 
Larmor  radius  in  the  poloidal  magnetic  field,  then  the  destabilizing 
effect  of  the  bootstrap  current  dominates  that  of  ion  inertia  for 
all  m.   Furthermore,  the  ion  acoustic  wave  (od  -  kn  c  )  is  also 
driven  unstable  by  the  bootstrap  current  in  this  case. 

A.  quasilinear  analysis  of  the  trapped  particle  diffusion 
resulting  from  these  instabilities  shows  that,  for  a  given  wave 

energy  spectrum,  the  particle  flux  (averaged  over  a  magnetic  surface) 

2 
is  proportional  to  I/Bq- 
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Neoclassical  Diffusion  in  Axisyrametric  Torus 
in  the  "intermediate  Region" 
Dilip  K.  Bhadra,  Chuan  S.  Liu,  and  Tihiro  Ohkawa 
Gulf  General  Atomic  Company 
San  Diego,  California 

Higher-order  collisional  effects  on  the  neoclassical  diffusion 
in  the  "intermediate"  region  of  a  toroidal  plasma  are  studied.   Since 
the  major  contribution  to  the  diffusion  is  from  the  Iow-Vm  particles, 

the  Fokker- Planck  equation  is  solved  for  an  asymptotic  solution  with 

1/2 
small  parameter  (v  /v, )  '      ( v  ,  v ,  being  the  collision  and  transit 

frequency  respectively).   The  resulting  diffusion  coefficient  is  of 

the  form  D  =  D  [l+a(v  /v   )-^^    -  ^{v   /v  ) . . .  ]  where 

O  C    L  C   Tv 

1/2  2 

D  =  (7r/2)  •  (2Trp  /lR)(p  Q,   )    is  the  result  previously  obtained  by 

Galeev  and  Sagdeev.   The  coefficient  a  is  small  and  vanishes  in  the 

absence  of  a  radial  electric  field.   The  coefficient  P  is  positive 

definite,  of  the  order  of  unity.   Thus  the  major  higher-order 

collisional  effects  tend  to  reduce  the  diffusion  coefficient  from 

D  .   This  decrease  is  due  to  the  reduction  of  the  random  walk 
o 

frequency.   The  half-pov/er  of  (v  /v.)  is  inherent  to  this  kind  of 
boundary- layer  type  problem. 


Work  supported  by  the  U.  S.  Atomic  Energy  Commission,  Contract 
No.  AT(04-3)-l67,  Project  Agreement  No.  38. 
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Shocks  and  Diffusion  in  Low-p  Toroidal  Plasmas 
N.  K,  Winsor,  E.  C.  Bowers,   M.  A.  Hellberg,    and  J.  M.  Dawson 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  N.J.  0854o 

The  low-P  fluid  equations  have  been  examined  for  axisyrametric 
solutions  which  are  stationary  (  ^^  =  0)  up  to  the  order  at  which 
resistive  diffusion  causes  the  macroscopic  parameters  to  change. 
With  an  ordering  of  the  resistivity  motivated  by  experimental 
parameters,  the  stationary  solutions  display  poloidal  rotation  and 
a  density  profile  with  a  discontinuity  curve (shock)  extending 
radially  from  the  magnetic  axis  to  the  wall.   When  viscosity  is 
added  to  the  equations,  both  shocked  solutions  and  unshocked 
solutions  are  possible,  depending  on  the  relative  magnitude  of  the 
resistivity  and  viscosity.   Numerical  simulation  studies  agree 
with  these  analytic  solutions.   Cases  with  more  viscosity  and  less 
resistivity  reveal  no  shock,  a  material  flux  approximately  classical 
(Pflrsch-Schluter)  and  a  negative  potential  on  axis  as  in  ST 
tokamak  experiments.   When  resistivity  is  lowered,  a  density  shock 
yields  fluxes  more  than  an  order  of  magnitude  greater  than 
classical  and  the  axis  is  positive,  as  in  Model  C  Stellarator 
measurements.   Movies  of  the  density  contours  for  these  cases 
illustrate  the  physical  processes  involved. 

Present  address:   Imperial  College,  London,  England 
On  leave  from  University  of  Natal.   The  author  gratefully 
acknowledges  receipt  of  an  Overseas  Study  Grant  from  the  South 
African  Atomic  Energy  Board. 
This  work  was  supported  by  U.S.  Atomic  Energy  Commission  Contract 
AT( 30-1) -1238;  Use  was  made  of  computer  facilities  supported  In 
part  by  National  Science  Foundation  Grant  NSF-GP  579- 
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stabilization  of  Dissipative  Trapped  Particle  Instability- 
David  W.  Ross  and  Marshall  N.  Rosenbluth 
The  Institute  for  Advanced  Study,  Princeton 

The  calculation  of  the  stability  of  the  trapped  particle 
drift  mode,  which  is  driven  by  electron  collisions,  is  amended 
to  include  untrapped  ion  Landau  damping  in  addition  to  the 

trapped  ion  collisional  damping  considered  previously.    It  is 

2 
found  that  the  Landau  term  is  stabilizing  if  dT/ir  =  0,   but 

that  it  changes  sign  if  d  In  T/d  In  r  >  2/3.   In  the  latter  case 

the  nixmerical  results  predict  that  ion  collisions  are  insufficient 

to  stabilize  the  mode  under  typical  operating  conditions  of 

several  proposed  Tokamaks. 

David  W.  Ross,  Marshall  N.  Rosenbluth,  and  D.  P.  Kostomarov, 
Bull.  Am.  Phys.  Soc .  15,  l400  (1970). 
^R.  Z.  Sagdeev  and  A.  A.  Galeev,  Sov.  Phys.-JETP  24,  1172(1967). 
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Computation  of  Finite-Beta  Plasma  Equilibria 
In  Minimum- B  Mirror  Systems 
David  V.  Anderson  and  John  Killeen 
University  of  California,  Lawrence  Radiation  Laboratory 

Plasma  equilibria  have  been  computed  for  various  containment 
schemes  which  share  the  common  feature  of  effective  two-dimensional 
symmetry.   We  present  a  method  for  computing  the  equilibria  in  the 
general  three-dimensional  minimum-B  mirror  system,  and  give  numerical 
results  for  finite-beta  equilibria  in  the  Baseball  II  magnetic 
configuration.   We  are  also  applying  this  code  to  the  2X  II 
experiment  and  to  the  design  of  the  mirror  fusion  experiment  (MFX) . 

For  the  case  of  open  field  lines  the  hydromagnetic  equilibria 

2 
equations  with  tensor  pressure  reduce  to  the  equation  V  9  =  l/v(V0  Vv) 

— *•  2 

where  B  =  l/v(Ve)and  v  =  (pn-Pj_)/B  -1.   We  can  write  9   as  ^^+^^' 

where  9      is  the  potential  function  of  the  vacuum  magnetic  field, 

i.e.,  B  -  -V9  ,  and  9      is  the  potential  function  of  the  magnetic 
c      c       p 

field  due  to  the  plasma.   We  solve  the  equation  for  9  (r, 9, z)  on  a 

cylindrical  domain  by  finite  difference  methods.   An  implicit 

iterative  algorithm  similar  to  the  ADI  scheme  Is  used  to  solve  the 

finite  difference  form  of  the  equation.   We  assume  that  the 

distortion  of  the  original  vacuum  field  by  the  plasma  is  negligible 

at  sufficient  distance  from  the  plasma  region.   This  distance 

determines  r    and  z    of  the  cylindrical  domain  and  on  this  outer 
maix      max         "^ 

boundary  9     --  0.   By  employing  a  sparse  mesh  in  the  exterior  region 

IT 

we  can  extend  the  boundary  surfaces  to  verify  the  effect  of  this 

condition.   The  location  of  the  coils  poses  no  problsm-the  boundary 

surface  can  enclose  the  coils-   We  use  the  pressure  functions  P|_(B) 

and  P||(B)  given  by  Taylor  and  obtain  solutions  corresponding  to  a 

sequence  of  betas  until  stability  is  violated. 

*Work  performed  under  the  auspices  of  the  U.S..  Atomic  Energy  Commission. 
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Limiter  Induced  Perturbations  and  the  Small  Gyro  Radius  Limit 
Donald  Dobrott  and  Harold  Grad 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

Although  Tokomaks  are  designed  to  have  approximately 
symmetric  fields,  the  limiter  (and  fringing  electric  fields) 
makes  the  geometry  strongly  asymmetric.   One  effect  is  to 
preferentially  scrape  off  circulating  as  compared  to  trapped 
particles.   More  significantly,  the  limiter  partially  blocks 
the  flow  of  toroidal  current.   The  simplest  non-dissipative 
macroscopic  model  which  allows  current  and  flow  to  depart 
from  magnetic  flux  surfaces  is  one  with  Hall  effect.   We 
examine  the  fully  non-linear  finite  P  theory  of  this  model 
in  axially  symmetry.   The  system  of  equations  is  eighth  order: 
twice  elliptic,  for  subsonic  flow  (distinct  boundary  conditions 
for  magnetic  and  fluid  flow  stream  functions),  plus  four  real 
characteristics  (flux  lines  and  streamlines,  each  counted  twice) 
The  limit  of  small  Hall  effect  (small  gyro- radius)  is  very 
singular  and  includes  unusual  boundary  layer  effects. 
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Low  Frequency  Drift  Cyclotron  Instability 
Burton  D.  Fried  and  Charles  F.  Kennel 
University  of  California,  Los  Angeles 

We  consider  the  instability  associated  with  ions  streaming 
across  a  magnetic  field  at  sufficiently  high  velocity,  V^,  that 
the  ion  cyclotron  radius  exceeds  the  plasma  dimensions. 

Forslund,  et  al  have  called  attention  to  the  possible 
importance  of  this  instability  for  plasma  heating,  and  have 
presented  solutions  of  the  linear  dispersion  relation  and  also 
the  results  of  particle  simulation  calculations  of  the  nonlinear 
behavior  for  the  case  of  purely  perpendicular  propagation  of  the 
unstable  waves  involved,  the  frequencies  being  of  order  electron 
cyclotron  frequency  (O  ) .   Recent  attempts  to  observe  this  in- 
stability in  laboratory  plasmas  were  carried  out  by  J.  M.  Sellen 
at  TRW  Systems  and  by  P.  Barrett  and  R.  Taylor  at  UCLA;  in  both 
cases,  it  was  found  that  the  dominant  instability  appears  to  be  of 

order  V„/r   =  (V_/a  )  O  <<  O  ,  due  to  the  fact  that,  as  in  all 
D  ce    ^  D  e'   e     e  ' 

laboratory  configurations,  k-, -,  is  not  exactly  zero.   (Typical 

minimum  values  of  k., .,  r   are  .01  to  .05).   We  show  that  the 

11  ce  ' 

instability  in  this  case  is  essentially  an  ion  acoustic  wave,  the 

dispersion  relation  being  modified  in  certain  essential  ways  by 

the  magnetic  field.   The  fastest  growing  waves  have  k.r   of  order 

1  and  7   =  Imw  of  order  (O  /O.)^^'^ . 

If  T  /T.  is  large,  then  the  behavior  is  governed  entirely  by 

the  two  dimensionless  parameters,  Q  =  (V-p^/a  )  (  V^k  r  )  and 

A  =  (O  /co   )  ,  where  a  =  /21   /m.   As  T  /T.  decreases,  the  growth 
^  e  pe'  '        e    ^   e         e   i  '  ^ 

rate  is  diminished  by  ion  Landau  damping.   Illustrative  curves 
showing  the  dependence  on  Q  and  A  and  comparison  with  the  experi- 
mental results  will  be  presented. 
■'"Forslund,  Morse  and  Nielson,  Phys .  Rev.  Lett.  25  1266  (1970) 


■X- 

Electron  Cyclotron  Drift  Instability 
D.  W.  Forslund,  R.  L.  Morse  and  C.  W.  Nielson 

Further  studies  of  the  electron  cyclotron  drift  instability 

indicate  that  it  is  this  instability  and  not  an  essentially  magnetic 

field  free  ion  acoustic  instability  which  is  responsible  for 

anomalous  diffusion  (or  resistivity)  in  high  density  containment 

experiments  and  in  collisionless  shocks.   This  instability  is  much 

less  sensitive  to  T  /T.  than  the  ion  acoustic  instability  and  in  fact 

persists  when  T.  >  T  and  V  «  V  (V,  and  V  are  the  drift  and 

electron  thermal  velocities  respectively) .   A  rough  instability 

criterion  is  V,/V   >  cd  /oO  and  V,  >•  V.  although  some  instability 
d  e    c   p      d  '^  1 

exists  for  smaller  drifts.   Detailed  od  vs  k  root  diagrams  will  be 
shown  as  well  as  namerical  simulations  of  the  resulting  electron 
cyclotron  turbulence.   In  particular,  anomalous  diffusion  will  be 
clearly  demonstrated.   The  implications  for  high  P  controlled 
fusion  devices  will  be  discussed. 


This  work  performed  under  the  auspices  of  the  U.  S.  Atomic 
Energy  Commission. 

""■D.  W.  Forslund,  R.  L.  Morse  and  C.  W.  Nielson,  Phys.  Rev. 
Letters,  25,  1270  (1970). 


-20- 


Nonlinear  Theory  for  Plasma  Driven  by  a  Large-Amplitude, 

High- Frequency  Field 
E.  J.  Valeo,  C.  Oberman,  and  ¥.  L.  Kruer 

By  means  of  an  asymptotic  multiple- time- scale  technique,  we  . 
systematically  develop  equations  capable  of  describing  the  non- 
linear evolution  of  a  parametrically  driven  homogeneous  plasma. 
The  small  expansion  parameter  is  the  electron-to-ion  mass  ratio. 
A  general  expression  for  the  nonlinear  (in  both  the  pump  and 
fluctuation  intensities),  high-frequency  energy  absorption  rate  is 
derived  as  a  function  of  the  fluctuation  amplitudes  which  may  be 
computed  nonlinearly  within  the  framework  of  the  theory.   The 
linearized  form  of  these  equations  is  compared  with  the  results 
of  others.   A  novel  consequence  of  the  linearized  theory  is  the 
prediction  of  a  finite  phase  relation  between  the  self-consistent 
high-frequency  waves  and  the  pump  wave,  which  necessitates  a 
modification  of  the  random-phase  approximation  of  usual  weak 
turbulence  theory.   Upon  consideration  of  the  evolution  of  the 
electron  distribution  function,  but  with  neglect  of  the  competitive 
processes  of  mode  coupling,  a  saturation  value  for  the  wave  amplitudes 
is  calculated  and  is  seen  to  be  in  reasonable  agreement  with  results 
from  a  numerically  simulated  plasma.   A  subsidiary  wave  vector 
expansion  of  the  equations  is  carried  out  and  the  important  processes 
of  electron-by-ion  wave  scattering  and  electron-to-ion  wave  decay 
are  shown  to  be  the  dominant  nonlinear  mode  coupling  processes. 
Upon  transformation  to  normal  mode  (Bloch  function)  form,  the 
equations  for  the  electric  potential  are  shown  to  describe  a  regime 
where  growth  and  oscillation  are  of  competitive  size.   The  results 
of  the  analysis  of  the  effects  of  these  important  mode-coupling 
processes  on  saturation  will  be  given. 

This  work  was  supported  by  U.S.  Atomic  Energy  Commission  Contract 
AT( 30-1) -1238;  one  of  the  authors,  E.J.V. ,  is  a  N.S.P.  predoctoral 
fellow.  -21- 


Suppression  of  the  Whistler  Instability  by 
Relativistic  Effects 

C.  L.  Hedrick,  Jr. 
Oak  Ridge  National  Laboratory,  Oak  Ridge,  Tennessee  37830 

The  anisotropy  driven  whistler  instability  has  been 
predicted  to  be  particularly  virulent  for  ECH  plasmas. 
Indeed,  it  has  been  observed  in  a  number  of  experiments.   In 
this  paper  it  is  shown  that  relativistic  effects  can  reduce  the 
growth  rate  to  a  negligible  level  as  a  consequence  of  the  relati- 
vistic spread  in  the  cyclotron  frequency,  which  causes  resonant 
particles  in  the  tail  of  the  distribution  function  to  contribute  to 
damping  rather  than  growth. 

For  several  distribution  functions  it  is  found  that,  for 
fixed  anisotropy,  the  growth  rate  decreases  dramatically  with 
temperature  (in  many  cases  the  growth  rate  becomes  negligible  for 
temperatures  which  are  small  compared  to  the  rest  energy  of  an 
electron) .   In  general,  distribution  functions  with  higher  tails 
show  a  more  pronounced  decrease  in  growth  rate.   These  results 
are  in  keeping  with  the  experimental  observations  of  N.  Lazar 
that  two- frequency  heating  simultaneously  Increases  the  tall  of 
the  distribution  function  and  suppresses  the  whistler  Instability. 
The  results  are  also  consistent  with  the  observations  of  R.  Dandl 
(no  whistler  instability  is  observed  and  the  distribution  function 
always  has  a  rather  high  tail) . 

Research  sponsored  by  the  U.  S.  Atomic  Energy  Commission  under 
contract  with  the  Union  Carbide  Corporation. 
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High  Frequency  A.  C.  Electrostatic  Plasma  Instabilities 
J.  P.  Freidberg  and  B.  M.  Harder 
Los  Alamos  Scientific  Laboratory 

The  instabilities  induced  in  a  plasma  under  the  influence 
of  an  electric  field  oscillating  at  near  the  plasma  frequency 
are  studied.   By  numerically  solving  the  integral  equation  which 
describes  the  linearized  behavior  of  the  plasma  we  are  able  to 
formulate  a  rather  detailed  and  comprehensive  picture  of  the 
stability  properties  for  a  wide  range  of  parameters. 

This  work  performed  under  the  auspices  of  the  U.  S.  Atomic 
Energy  Commission 
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Dynamic  Stabilization  of  Drift  Waves  by  Radial  and 

-X- 

Azimuthal  High- Frequency  Fields 
P.  K.  Kaw  and  Y.  C.  Lee 
Plasma  Physics  Laboratory,  Princeton  University, 
Princeton,  N.  J.  085^0 

We  have  investigated  the  stabilization  of  density  gradient 

driven  low-frequency  drift  instability  (w  <<  ^^;    k^^a^jk^^a^  <<  1; 

a  and  a.  refer  to  the  electron  and  ion  Larmor  radii  and  the 
e      1 

rest  of  symbols  have  usual  meanings)  and  the  drift  cyclotron 
instability  (<*)  ~  ^.,    k  a.  ~  1)  by  radial  or  azimuthal  high- 
frequency  fields  with  a  frequency  of  the  order  of  upper  hybrid 
resonant  frequency.   This  scheme  of  stabilization  is  much  more 
efficient  than  the  conventional  method  of  applying  axial  high- 
frequency  fields  because  the  drift  waves,  which  essentially 
propagate  at  right  angles  to  the  magnetic  field,  couple  more 
strongly  to  fields  oriented  in  that  direction.   Furthermore,  one 
can  also  stabilize  the  more  dangerous  flute-like  drift  cyclotron 
waves  (with  kn  =  0)  by  this  scheme  in  contrast  to  the  conventional 
method  using  axial  high-frequency  fields  only. 


-""Ya.  B.  Fainberg  and  V.  D.  Shapiro,  Sov.  Phys .  JETP  25,  189  (1967); 
M.  Okamoto,  T.  Amano,  and  K.  Kitao,  J.  Phys.  Soc.  Japan  29,  104l(1970) 
This  work  was  supported  by  U.  S.  Atomic  Energy  Commission 
Contract  AT(30-1)-1238. 
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High  Pressure  Tokamak  EquilibriTjm 

H.  R.  Strauss 

University  of  Texas  at  Austin 

Assuming  an  ideal  MHD  equilibrium,  we  find  there  is  no 

limit  to  the  maximum  pressure  p    obtainable  in  a  Tokamak  of 

^  max 

given  magnetic  field  strength  B,  aspect  ratio  e   ,  and  rotational 

transform  i.   There  is  a  limit  to  the  total  energy  a  Tokamak 

-    -  2  2 

can  contain;  this  may  be  expressed  as  p  <  p  =  ae{  i/2{  £/2Tr)    B  /Sir 

where  p  is  the  pressure  averaged  over  the  plasma  volume  and  a  is  a 
number  of  order  unity.   The  pressure  and  rotational  transform 
have  their  maximum  values  on  the  magnetic  axis.   As  p  and  p^^^ 
increase  with  £,    e,    B  fixed,  the  distance  B  from  the  magnetic 
axis  to  the  outer  edge  of  the  plasma  decreases.   When  5  <  .Jr, 
where  r  is  the  plasma  radius,  the  toroidal  current  begins  to 
reverse  on  the  inner  edge  of  the  plasma.   For  6  <  .6r,  p  has 
reached  its  limiting  value  p  ,  and  &  =  (3p  /Pmax^^'  '^^^   axis 
moves  to  the  edge  in  the  limit  p   -♦  oo .   In  general,  there  are 
two  classes  of  high  p  Tokamak  equilibria,  "deep"  solutions  with 

arbitrarily  large  p    and  "shallow"  ones  with  p   ,  limited  to 
•^  ^     ^max  ^msLX 

about  2p  .   The  magnetic  axis  of  shallow  solutions  does  not 

displace  very  much  as  p  and  p    tend  to  their  limiting  values. 

max 

In  deriving  all  these  results,  p(VO  and  I  {f)    have  been  chosen 
quadratic  in  f,    where  p  is  the  pressure,  I  is  the  poloidal  current 
flux,  and  t   is  the  poloidal  magnetic  flux.  (For  p(7/')  linear  in  f, 
there  are  only  shallow  solutions) . 
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Thermonuclear  Power  Production  by  Nonmaxwellian 
Ions  in  a  Closed  Magnetic  Field  Configuration 
J.  M.  Dawson,  H.  P.  Furth,  and  F.  H.  Tenney 
Plasma  Physics  Laboratory,  Princeton  University 
Princeton,  New  Jersey  0854o 

High-energy  neutral  beam  injection  into  the  tokamak  configu- 
ration can  be  used  to  implement  the  ancient  idea  of  a  target- 
plasma  fusion  reactor.   In  a  cold-triton  plasma  with  T  =  5  keV, 
a  l80-keV  deuteron  on  the  average  2.4-folds  its  energy  while 
thermalizing.   (This  factor  includes  initial  kinetic  energy, 
fusion  release,  and  standard  blanket  reactions.)   Assuming 
moderate  efficiencies  of  conversion  from  thermal  to  electrical 
power,  net  power  can  be  produced.   Optimization  calculations 
have  been  carried  out  and  will  be  presented.   The  main  advantage 
of  the  scheme  is  to  lower  the  appropriate  Lawson-type  requiremnets 
significantly  relative  to  the  one-component  reactor  case.   The 
stability  problem  will  also  be  discussed:   it  does  not  appear  to 
be  severe,  if  the  energetic-ion  component  is  fairly  isotropic. 


* 
This  work  was  supported  by  U.S.  Atomic  Energy  Commission  Contract 


AT(30-1)-1238. 
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Thermonuclear  Reaction  Waves  at  High  Densities 
Ming-sheng  Chu 
Colimahia  University 

Phenomena  are  studied  of  the  ignition  and  the  development 
of  plane  thermonuclear  reaction  waves  propagating  through  solid 
density  deuterium-tritium  plasmas.   The  state  of  the  plasma  is 
described  by  a  two  temperature  gas  and  optically  thin.   The 
equations  governing  the  change  of  density,  velocity,  electron 
temperature,  ion  temperature  and  degree  of  reaction  are  derived 
from  the  conservation  laws  and  the  laws  of  reaction  kinetics. 
The  physics  are  contained  within  a  set  of  five  coupled  non-linear 
partial  differential  equations.   Terms  in  the  equations  take 
explicit  account  of  thermonuclear  reactions  and  heating,  conduction, 
electron-ion  equilibration,  bremsstrahlung  cooling  and  other  fluid 
dynamic  effects.   Qualitative  behavior  of  the  physical  variables  is 
discussed.   Solutions  of  the  equations  are  obtained  by  employing 
finite  difference  techniques.   Due  to  the  properties  of  the  plasma 
and  the  form  of  the  reaction  cross  section,  substantial  physical 
differences  exist  between  the  fusion  driven  reaction  waves  and 
detonation  waves  in  chemically  reactive  gases.   The  ignition  energy 
threshold  is  obtained  for  simple  geometries  and  the  above  threshold 
case  is  studied.   Throughout  the  development,  the  wave  structure 
changes  from  pure  conduction,  to  plasma-schock-like  behavior,  and 
finally,  it  is  dominated  by  thermonuclear  heating  and  bremstrahlung 
cooling  effects.   The  final  steady  wave  structure  is  presented. 
The  thermonuclear  energy  yield  from  small  D-T  pellets  has  also 
been  studied. 

*¥ork  supported  by  AFOSR  contract  AF44620-71-C-0100 
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Multi-Species  Fokker  Planck  Calculations  and  Q  Values 
For  D-T  and  D-?Ie-^  Mirror  Reactors 
J.  P.  Holdren,  A.  H.  Futch,  Jr.,  J.  Killeen,  A.  A,  Mirin  and  R.  F.  Post 
University  of  California,  Lawrence  Radiation  Laboratory 

Livermore,  California 

The  feasibility  of  the  mirror  machine  as  a  fusion  reactor 
depends  in  large  part  on  the  achievable  values  of  Q,  the  ratio  of 
thermonuclear  power  to  injected  power  in  steady  state.   Under  the 
assiimption  that  the  end  losses  are  due  primarily  to  weak  Coulomb 
collisions,  Q  can  be  determined  from  a  solution  of  the  coupled, 
time-dependent,  Fokker- Planck  equations  for  the  species  present. 
The  computer  code  used  in  the  present  work  obtains  distributions 
f.(v,  t),  for  an  arbitrary  nixmber  of  species;  the  mirror  machine  is 
approximated  by  a  magnetic  square-well,  and  the  *  and  9   dependences 
of  the  distributions  in  velocity  space  are  removed  ty  assuming 
azimuthal  symmetry  and  by  separation  of  variables,  respectively. 
Improvements  over  previous,  related  work  are:   explicit  treatment  of 
3-component  systems  such  as  D-T-e  and  D-He  -e;  refined  calculation  of 
reaction  rates,  <"  ,  v>,  integrating  explicitly  over  numerically 
determined  f .  for  both  reactants;  and  more  rigorous  treatment  of 
the  Coulomb  logarithm,  InA.   Results  to  be  presented  emphasize 
Q  values  for  the  D-He-^  fuel  cycle,  with  injection  energies  in  the 
ranges,  200  keV  j^  Ej^  j<  800  keV,  800  keV  _<  Ej^^  <   2000  keV.   Values 
of  Qy,,p  versus  injection  energy,  their  sensitivity  to  the  means  of 
calculating  «tv>,  and  new  results  concerning  the  mass  dependence 
of  individual  nx ' s  will  also  be  discussed. 


■X- 

Work  performed  under  the  auspices  of  the  U.  S.  Atomic  Energy 
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¥111  Ion  Cyclotron  Waves  Heat  Tokamak  Plasmas? 

J.  M.  Kindel  and  F.  W.  Perkins 
Plasma  Physics  Laboratory,  Princeton  University, 
Princeton,  N.  J.  0854o 

The  object  of  this  work  is  to  determine  the  conditions  under 
which  ion  cyclotron  waves  can  efficietnly  heat  ions  in  the  central 
regions  of  tokamak  plasmas,  and  what  consequences  (such  as  electron 
heating)  ensue  if  the  ion  cyclotron  waves  are  not  launched  in  the 
optimum  manner.   Our  model  of  a  tokamak  is  a  plasma  slab  of  varying 
density  perpendicular  to  a  confining  magnetic  field  with  a  varying 
magnitude  but  constant  direction.   As  a  result,  co  and  kn  remain 
fixed  as  a  wave  propagates  and  the  dispersion  relation  determines 

the  (complex)  wave  vector  k  .   Typical  results  are  that  ion  cyclotron 

2       2      l/^         —1 
waves  should  be  launched  with  kn  x    [ca   .    Q,  /c   v.,  .  ]  '^ -*  a  0.7  cm   and 

11     pi  c    thi 

a   «  0.95(^  where  Q,     is  the  ion  cyclotron  frequency  at  the  center  of 
the  plasma.   If  kn  is  too  large,  the  ion  cyclotron  waves  are  strongly 
evanescent  and  do  not  reach  the  central  plasma.   For  values  of  kn 
somewhat  too  small,  the  ion  cyclotron  wave  converts  into  the  slow 
Alfven  wave  which  principally  heats  electrons.   For  very  small 
values  of  kn,  there  is  a  mode  conversion  into  an  electrostatic  wave 
at  the  plasma  surface  with  possible  surface  heating. 


This  work  was  supported  by  U.  S.  Atomic  Energy  Commission 


Contract  AT(30-1)-1238. 
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Cross  Field  Injection  of  a  Relativistic  Electron  Beam 
into  a  Magnetized  Plasma 
Roswell  Lee  and  R.  N.  Sudsin 
Cornell  University  and  Naval  Research  Laboratory 

In  the  Astron  experiment  and  in  a  similar  experiment  at 

Cornell  University  ,  relativistic  electrons  are  injected  into  a 

magnetic  mirror  to  form  a  layer  of  circulating  electrons.   When 

the  beam  of  electrons  is  injected  into  a  preionized  plasma  it  is 

magnetically  neutralized  by  a  return  current  which  is  induced  in 

the  plasma  by  the  injection  process.   We  have  computed  the  magnitude 

and  decay  time  of  this  return  current  on  the  assumption  that  the 

plasma  is  cold  and  the  beam  is  rigid.   We  find  that  such  a  return 

will  flow  mostly  within  the  beam  cross-section  when 

cxo)  (1 +cd~/(jd'")  '    /C  >>  1,  where  a  is  the  beam  radius,  co  is  the 
p^    c  p'    '      '  p 

plasma  frequency,  and  'jd  is  the  cold  electron  cyclotron  frequency. 

The  return  current  flows  for  a  time  of  the  order  of 

0)  .    [u,   ocD  /v  71,  where  p.  =  m  /m.,  v  is  the  beam  velocity, 
ci   "-      p  o  e  i'   o  ^ 

2   ^  —1/2 
7  =  (1-v  /c~)    ,  and  cd  .  is  the  cold  ion  cyclotron  frequency. 


H.  H.  Fleischmann,  H.  Davitian,  R.  E.  Kribel,  and  J.  A.  Nation, 
Bull.  Am.  Phy.  Soc.  15,  1^59- 
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Precession  of  a  Strong  E-Layer 
R.  V.  Lovelace  and  R.  N.  Sudan 
Cornell  University  and  Naval  Research  Laboratory 

Observations  of  weak  relativistic  E-layers  in  the  Astron 
indicate  that  the  motion  is  rigid  in  accord  with  the  "super 
particle"  model  of  Purth  (1)  in  which  electrons  of  the  layer  are 
assumed  to  precess  as  a  rigid  body  with  a  frequency  equal  to  the 
precession  frequency  of  a  single  electron  in  the  external  fields. 
By  perturbation  of  the  Vlasov  equation  for  an  E-layer  of  arbitrary 

strength  we  investigate  the  affect  the  self-magnetic  field  has  on 

2  /  2 
the  precessional  motion.   In  terms  of  the  parameter  I  =  co  /2fi 

(where  -iJ-p,  is '  the  relativistic  electron  plasma  frequency  and  Q, 

the  cyclotron  frequency)  we  find  that  the  precession  frequency  is 

reduced  from  that  for  a  single  electron  by  a  factor  (1+1).   In  the 

weak  E-layer  limit,  |  <<  1,  the  precession  frequency  agrees  with 

that  of  Furth,  but  for  a  strong  E-layer  with  loading  factor 

C  ~  1,  I  approaches  r/2a,  where  r  is  the  layer  radius  and  a  its 

thickness.   Non-rigid  motion  of  the  layer  appears  in  our  analysis 

as  a  damping  of  the  precession  in  the  absence  of  external  forces 

on  the  layer;  the  damping  is  negligible  and  the  motion  is  rigid 

for  large  enough  self-fields,  i    >>  Ap/p  where  Ap/p  is  the  fractional 

spread  in  canonical  momentum.   Stability  of  the  precession  mode 

depends  on  the  sign  of  the  mode  energy,  the  positive  energy  case 

being  stable  under  the  influence  of  external  dissipative  forces; 

the  sign  of  the  mode  energy  is  shown  to  be  opposite  to  that  of  the 

external  field  index  independent  of  the  layer  strength.   For  a  weak 

E-layer  with  a  fraction  f  of  its  charge  unneutralized,  the  precession 

is  stable  only  for  (fC  +n)  <  0,  where  n  is  the  field  index. 


-^H.  P.  Purth,  Phys.  Fluids  8,  2020  (1965) 
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toroidal  field  effects  e.g.  trapped  particle  instabilities. 
The  relativistic  proton  E-layer  currents  substitute  for  the 
plasma  current  required  to  create  the  shear  thus  allowing  a 
continuous  operation  in  contrast  to  the  pulsed  operation  of 
the  Tokamak.  At  the  same  time  the  E-layer  particles  will 
heat  the  plasma  to  fusion  temperature. 
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E-Layer  Equilibria  in  the  Presence  of  a  Toroidal  Field 

N.  C.  Christofilos 
Lawrence  Radiation  Laboratory,  University  of  California 
Livennore,  California  9^550 

The  addition  of  toroidal  field  to  improve  the  plasma  P  in 
plasma  confined  in  the  Astron's  E-layer  was  suggested  several 
years  ago.    It  was  only  recently  however  that  the  effect  of  the 
toroidal  field  in  the  E-layer  equilibrium  has  been  investigated. 
The  investigation  so  far  is  limited  to  a  simplified  model  of  the 
E-layer  represented  by  a  straight  beam  of  circular  cross  section 
while  the  toroidal  field  is  parallel  to  the  beam  axis.   There  are 
two  classes  of  solutions  depending  whether  the  E-layer  is  paramagnetic 

or  dlamagnetic  to  the  toroidal  field.   The  first  class  of  solution  was 

p 
discovered  by  S.  Yoshikawa.    In  this  solution  the  relativistic 

electrons  follow  paths  almost  parallel  to  the  magnetic  field  thus 

minimizing  synchrotron  radiation  losses.   The  solution  given  in 

Reference  2  is  accurate  for  large  rates  of  v,  where  v  =  (l/l7,000  7), 

I   is  the  beam  current  and  7   is  the  electron  energy  expressed  in  rest 

energy  coils.   A  more  general  solution  valid  for  any  value  of  (v) 

will  be  presented. 

In  the  second  class  of  solution  the  E-layer  is  diaraagnetic  to 

the  externally  applied  toroidal  field.   In  this  way  the  stability 

properties  of  the  E-layer 's  closed  magnetic  well  are  combined  with 

Tokamak's  shear  to  provide  an  improved  plasma  confinement  hence 

higher  P.   In  addition  the  E-layer  currents  generate  a  depression 

or  a  minimum-B  in  the  toroidal  field,  thus  eliminating  all  the 


Work  performed  under  the  auspices  of  the  U.  S.  Atomic  Energy 

Commission. 

T.  K.  Fowler,  et  al.,  Conf.  on  Plasma  Physics  and  Controlled  Nuclear 

Fusion  Research,  Novosibirsk,  1968  (IAEA,  Vienna,  1,    981  (1969). 

2 
S.  Yoshikawa,  Physical  Review  Letters  26_,  295  (1971). 
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Formation  and  Interaction  of  Ion-Acoustic  Solitary 
Waves  in  a  Collisionless  Warm  Plasma 
P.  H.  Sakanaka,  T.  C.  Marshall,  and  C.  K.  Chu 
Plasma  Lab.,  Columbia  University,  New  York,  New  York  10027 

The  formation  and  interaction  of  solitary  waves  in  a  warm- 
ion,  hot  electron  plasma  have  been  studied  with  the  Vlasov 
equation  for  ions  and  Boltzmann  distribution  for  electrons-, 
with  the  purpose  of  determining  those  properties  hitherto  known 
only  for  cold  ions  and  the  K-dV  equation.   A  parallel  study 
has  also  been  made  using  the  fluid  equations  for  ions  with  a 
scalar  pressure  added. 

First,  by  using  the  warm  ion  fluid  equation,  with  isothermal 
electrons  given  by  Boltzmann  distribution,  we  have  shown  that 
for  the  same  amplitude, the  width  of  the  solitary  wave  is 
smaller  than  that  obtained  from  the  Korteweg-de  Vries  equation 
while  the  Mach  n-jmber  is  slightly  larger.   The  limiting  value 
of  the  amplitude  is  considerably  smaller  than  that  for  the 
case  of  the  cold  ions,  for  the  electron-ion  temperature  ratio 
T  =  30,  the  limiting  amplitude  is  only  1/3  of  that  for  the 
case  of  T  =  00  . 

Second,  with  ions  described  by  one-dimensional  Vlasov 
equation  and  electrons  by  the  Boltzmann  distribution,  we  have 
shown  that  an  arbitrary  propagating  pulse  in  a  plasma  always 
decays  into  solitary  waves.   We  have  studied  the  interactions 
of  two  solitary  waves,  a  large  wave  catching  up  with  a  smaller 
one.   We  have  found  that:   (1)  if  the  initial  amplitude  ratio 
of  the  two  waves  is  not  too  large,  the  two  waves  exchange  their 
amplitudes  and  bounce  apart;  (2)  if  the  initial  amplitude  is 
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large,  larger  wave  first  absorbs  the  smaller,  and  later  reemits 
it  in  the  back,  and  the  initial  amplitudes  are  conserved.   These 
properties  agree  qualitatively  with  the  Korteweg-de  Vries  re- 
sults, but  the  solution  of  the  ^Tlasov  equation  gives  a  more 
complete  picture  of  the  physics  of  the  interaction. 


Supported  by  Atomic  Energy  Commission  Contract  AT( 30-1) -395^ 
and  National  Science  Foundation  Grant  GK  1391X2. 


-35- 


Evolution  of  Shock  and  Nonlinear  Compression  Waves 

Pung  Nien  Hu 

Space  Sciences  Division,  Whittaker  Corporation 

Waltham,  Massachusetts  02154 

The  initial  value  and  boiindary  value  problem  of  weak  nonlinear 
plasma  waves  in  the  presence  of  a  perpendicular  magnetic  field  is 
studied,  based  on  a  system  of  moment  equations  of  Boltzmann 
equations  for  electrons  and  ions.   This  system,  rigorously  obtained 
by  proper  scalings  in  time  and  space  with  respect  to  the  wave  strength, 
describes  a  fully-ionized  plasma  with  all  possible  dissipative 
mechanisms  including  electron  and  ion  viscosity,  coupled  electron 
and  ion  heat  conduction,  thermal  diffusion  and  resistivity,  as  well 
as  all  the  distinct  Hall  effects.   In  all  the  plasma  parameter 
ranges,  it  is  found  that  the  evolution  of  weak  nonlinear  waves  is 
governed  by  the  Burgers  equation,  subject  to  given  initial  and 
boundary  conditions. 

The  formation  of  a  shock  and  the  evolution  of  a  compression 
pulse  are  studied  in  detail.   The  transit  time  and  distance  for  a 
fully  developed  shock  are  obtained.   In  some  parameter  ranges,  the 
transit  time  is  smaller  thain  the  electron  collision  time  while  the 
transit  distance  is  smaller  than  the  mean  free  path. 

The  time  scale  for  the  evolution  of  a  compression  pulse 
depends  strongly  on  its  initial  width.   A  shock-like  wavefront  will 
shortly  appear  when  the  initial  width  is  large  as  compared  with  the 
thickness  of  a  steady  shock.   (In  some  parameter  ranges,  this  shock 
thickness  is  smaller  than  the  mean  free  path.)   Eventually,  hovjever, 
a  compression  pulse  becomes  self-similar,  with  its  amplitude  decaying 
as  1/  \/r  and  its  width  spreading  as  v/T- 


-36- 


The  inclusion  of  the  electron  Inertia  effect  leads  to  the 
Korteweg-de  Vries-Burgers  equation  which  covers  a  wide  range 
of  wave  profiles. 


-X- 

Work  supported  by  the  Air  Force  Office  of  Scientific  Research. 
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Noise  Broadening  of  Trapped-Particle  Echoes 
E.  A.  Williams  and  C.  Oberman 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  N.J. 

It  has  recently  been  proposed  by  Liu  and  Wong  that  the 
behavior  of  magnetically  trapped  particles  can  conveniently  be 
studied  by  investigating  echoes  associated  with  them.   We  extend 
the  analysis  of  Liu  and  Wong  to  include  the  effect  of  interparticle 
collisions  or  weak  turbulence.   It  is  well  known  that  echo  phenomena 
in  general  are  sensitive  to  noise,  because  of  their  dependence  on 
the  preservation  of  fine  structure  in  the  distribution  function. 

Information  about  the  diffusion  of  magnetically  trapped 
particles  is  of  particular  interest  due  to  their  apparent  role  in 
the  stability  and  transport  of  toroidal  plasmas.   VJe  show  that  the 
effect  of  noise  on  the  echo  pulses  is  to  broaden  them,  provided 
the  effective  collision  time  is  long  compared  with  a  typical  bounce 
time  in  the  well.   If  this  condition  does  not  hold,  the  echo  is 
destroyed.   To  calculate  the  trapped-particle  distribution 
responsible  for  the  echoes,  we  require  the  propagator  for  a  particle 
trapped  in  a  magnetic  well.   This  we  compute  and  then  evaluate  the 
propagator  for  the  case  of  tokamak  geometry,  and  finally  specialize 
to  deeply  trapped  particles.   We  then  use  the  previous  results  to 
solve  the  echo  problem  and  discuss  the  results.   There  are 
interesting  magnetospheric  applications  as  well. 

* 
This  work  was  supported  jointly  by  U.  S.  Atomic  Energy  Commission 

Contract  AT( 30-1) -1238  and  U.  S.  Air  Force  Office  of  Scientific 

Research  Contract  P44620-70-C-0033. 
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Nonlinear  Interaction  of  a  Small  Cold  Beam  and 
a  Plasma,  Part  II  (Random  Phase  Approximation  Revisited) 
T.  M.  O'Nell  and  J.  H.  Winfrey 
University  of  California  at  San  Diego 
La  Jolla,  California 

In  a  previous  paper  ,  a  single  wave  model  is  used  to  describe 
the  initial  nonlinear  interaction  of  a  small  cold  beam  and  a  plasma. 
This  description  is  based  on  the  observation  that  after  several 
e-foldings  the  band  width  of  the  growing  spectrum  of  waves  is  so 
narrow  that  the  electrons  interact  with  a  very  nearly  pure  sinusoidal 
field.   However,  this  solution  is  only  locally  valid,  since  the 
small  remaining  spectral  width  produces  a  slow  spatial  modulation 
in  the  phase  and  amplitude  of  the  sinusoidal  wave.   This  spatial 
modulation  is  random,  since  the  phases  of  the  waves  within  the 
small  spatial  width  are  random. 

In  the  present  paper  a  nonlocal  solution  is  constructed  by 
giving  the  locally  valid  sinusoidal  solution  a  slow  spatial 
modulation  in  phase  and  amplitude.   Spatial  Fourier  transforms  oA'-er 
this  nonlocal  solution  are  evaluated  using  the  theory  of  random 
noise.   This  procedure  does  not  Involve  the  usual  assumptions 
associated  with  the  random  phase  approximation,  since  here  the 
dynamics  is  worked  out  completely,  before  any  averaging  is  done. 

"Nonlinear  Interaction  of  a  Small  Cold  Beam  and  a  Plasma", 

T.  M.  O'Nell,  J.  H.  Winfrey  and  J.  H.  Malmberg,  to  be  published  in 

Phys .  Fluids . 

Supported  in  part  by  the  National  Science  Foundation,  Grant  GP-27120, 
and  in  part  by  the  U.  S.  Atomic  Energy  Commission,  Contract  AT(o4-3)-3^ 
PA  85-13. 
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On  the  Existence  of  Force  Free  Magnetic  Fields  In  a  Simply 

Connected  Domain 
M.  Blneau 
Association  Euratom-Cea 
France 

A  force  free  magnetic  field  B(x)  carries  a  constant  current 

density  a(x)  along  each  magnetic  line.   This  fact  suggest  a  natural 

boundary  value  problem  for  a  force  free  field  In  a  simply  connected 

domain  D  with  boundary  Sd.   Let  S  be  some  part  of  ^D  and  D(B, S)  be 

the  part  of  D  connected  with  S  by  a  line  of  the  field  3.   The 

problem  Is  to  find  a  field  'b('x)  such  that  V.'b  =  0  In  D,  "B'Vo  =  0  and 

V  X  ^  =  PoB  In  D(B,  S),V  x  ^  =  0  In  D-  D(B,  S)  with  the  boundary  data 

of  B  on  3d  and  a  on  S.   (That  Is  n« B  =  B  given  where  n  Is  the  outer 

normal  of  ^D);  P  Is  a  parameter  that  measures  the  distortion  from  a 

vacuvun  field  "b  ,    solution  of  V«'b  =  0,  V  x  "b  =0,  'n-'s     ■=   B^  on  bD. 
o  o  o         o    n 

It  Is  assumed  that  P  <<  1.   With  j(x)  =  a(x)B(x)  in  D(B, S)  and 
j(x)  =0  in  D-D(B,  S)  the  solution  has  the  representation  B  =  B  +  PB  , 
^,  =  7^-  /v  i  X  "3"'  dr*  +-JP  /v  -  X  T'  ds'  (f '  means  f('x')  and  "x'  is 
the  variable  of  integration)  where  7   is  the  solution  of  the  Integral 
equation  on  3d  "7+-^  /n"  x  (V  i  x  7')dS'  =  "  "^  ^  Tv  i  x  "J*  dx ' . 
This  form  was  given  by  H.  Villat  to  the  solution  of  the  vortex 
problem  of  Poincare-Steklov  (in  terms  of  j  and  B  this  problem  amounts 
to  solve  the  same  equations  for  j(x)  given  in  D  and  B  =  0  on  3D)  . 
This  representation  allows  to  show  by  expansion  the  existence  of  a 
solution  when  the  field  B  is  regular  (does  not  vanish  in  D)  for 
small  enough  p.   If  B  has  unlformely  A-holder  continuous  first 
derivatives  and  a   has  unlformely  bounded  derivatives  on   3d  then  B 
has  unlformely  v-holder  continuous  first  derivatives  in  D  (v  <  A) . 

"""R.  Kress-Arch- Rational  Mech.  anal.  30  (1968)  33l-400. 
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Containment  in  Open-Ended  Systems  in  the  Presence  of 
Fluctuating  Electric  Fields 
R.  E.  Aamodt,  Department  of  Astro-Geophysics 
University  of  Colorado,  Boulder,  Colorado 

The  lack  of  constancy  of  the  energy  and  magnetic  moment  due  to 
fluctuating  electric  fields  need  not  drastically  effect  the  containment 
properties  of  mirror  traps.   If  the  amplitude  of  the  fluctuating 
electric  fields  are  small  enough  so  that  a  particle  E  x  B  drifts  only 
a  fraction  of  a  wavelength  in  a  cyclotron  period,  then  in  addition  to 
this  drift  the  main  effect  of  the  electric  fields  is  to  slowly 
modulate  the  Larmor  motion.   This  slow  modulation  can  "be  described 
by  the  equations  of  motion  averaged  over  a  Larmor  cycle.   In  the 
important  case  that  the  fields  have  (a)  frequency  close  to  a  harmonic 
of  the  gyro  frequency  (t>)  long  wavelength  along  the  external  magnetic 
field  and  (c)  short  wavelength  across  it,  these  time  averaged  equations 
are  shown  to  have  Hamiltonian  form  and  are  periodic.   Hence,  an 
adlabatic  invariant  for  the  perpendicular  motion  exists,  this  invariant 
reducing  to  the  usual  magnetic  moment  in  the  limit  when  the  fluctuating 
electric  fields  are  zero.   Motion  along  the  lines  of  force,  for  the 
averaged  equations,  implies  a  second  constant  of  the  motion  which 
corresponds  to  the  time  averaged  energy  in  the  frame  of  reference 
rotating  with  the  cyclotron  frequency.   Using  these  two  constraints, 
single  particle  ion  motion  in  a  mirror  system  is  examined  and  it  is 
shown  that  small  amplitude  waves  induce  large  magnetic  moment 
excursions.   This  implies  that  unstable  cyclotron  waves  which  depend 
sensitively  on  the  recurrent  phase  relations  of  the  simple  Larmor 
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cycle  will  be  stabilized  at  low  amplitudes.   However,  these 
constraints  also  suggest  that  when  the  fluctuating  fields 
vanish  well  inside  the  magnetic  mirror  field  maximum, 
particle  containment  in  the  presence  of  quasi-stationary  waves  is 
affected  only  as  an  ion  transits  a  "beach  region"  where  the 
frequency  of  the  waves  is  extremely  close  to  an  integral  multiple 
of  the  ion  cyclotron  frequency.   This  is  a  result  of  the  fact 
that  only  in  these  resonance  regions  is  the  particle  motion 
"irreversible."  Such  a  local  process  minimizes  the  overall  effect 
of  the  large  magnetic  moment  excursions  and  greatly  reduces  expected 
particle  loss  rates. 

Work  supported  in  part  by  the  U.  S.  Atomic  Energy  Commission 
under  Contract  AT(ll-l)-203l. 
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Effect  of  Field  Asymmetries  on  Neoclassical  Confinement 

in  Tokamaks 
M.  N.  Rosenbluth,  P.  H.  Rutherford,  and  R.  Hazeltine 
Plasma  Physics  Laboratory  and  the  Institute  for  Advanced  Study 

Princeton,  New  Jersey  085^0 

Experiments  are  being  carried  out  on  the  Princeton  ST  tokamak 
on  the  influence  on  confinement  of  local  increases  in  toroidal  field. 
¥e  have  considered  the  effect  of  such  a  symmetry  on  neoclassical  trans- 
port theory.   The  usual  axisymmetric  field  is  used  but  a  local  increase 
in  field  8B  is  simulated  by  imposing  specular  reflection  at  <t>  =  0 
for  particles  with  e  -  |-i.B„  <  i-L6B;  here  e  and  m-  are  the  energy  and 
magnetic  moment,  and  B  =  B„[l-  (r/R)  cos  6],      In  contrast  to  the  axi- 
symmetric problem  it  is  found  necessary  to  solve  for  a  curious  boundary 
layer  between  trapped  and  untrapped  particles:   this  solution  is 
developed  by  the  Wiener-Hopf  method.   In  the  low  collision  frequency 

regime  both  the  radial  diffusion  and  the  transport  due  to  the  driving 

2 
E- field   (trapped  particle  pinch  effect)  are  strongly  modified. 

Surprisingly,  in  the  case  of  small  rotational  transform  and 

8B/B  >>  r/R,  the  latter  becomes  a  particle  flux  of  order 

F  ~  (nE/B)(r/R)(B/8B)^'^^  (v^vR)-"-/^,  being  therefore  of  increased 

magnitude  at  low  collision  frequency,  and  radially  outward  in  direction. 


W.  Stodiek,  K.  Bol,  H.  Eubank,  and  S.  von  Goeler  (to  be  published) 

2 
P.  H.  Rutherford,  L.  M.  Kovrizhnikh,  M.  N.  Rosenbluth,  and 

F.  L.  Hinton,  Phys .  Rev.  Lett.  25,  IO9O  (1970). 
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This  work  was  supported  by  U.  S.  Atomic  Energy  Commission 
Contract  AT(30-1)-1238. 
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A  Lagrangian  Theory  for  Nonlinear  Wavepackets  in  a  Collisionless  Plasma 

R.  L.  Dewar 
University  of  Maryland 

A  rather  general  method  for  treating  the  nonlinear  interaction 
between  a  wavepacket  and  a  collisionless  plasma  is  presented. 
The  method  is  applicable  to  a  v;ide  class  of  waves,  including  laser 
pulses,  but  the  specific  wave  investigated  is  the  electron  plasma 
wave.   The  method  consists  in  representing  the  potential  and 
particle  disturbances  in  terms  of  slowly- varying  (in  space  and  time) 
complex  amplitudes  for  the  fundamental  and  harmonics,  substituting 
these  into  the  Low  Langrangian  for  the  plasma,  averaging,  and 
simplifying  by  elimination  of  all  amplitudes  but  that  of  the 
fundamental  of  the  wave  potential.   The  resultant  Lagrangian  gives 
rise  to  a  nonlinear  Schrodinger  equation  for  the  wave  potential. 
The  nonlinearity  arises  both  from  the  purely  nonlinear  frequency 
shift  and  from  quasilinear  modification  of  the  linear  dispersion 
relation.   The  latter  effect  includes  nonlinear  Landau  damping, 
which  is  shown  to  lead  to  instability  of  the  wave  envelope,  thus 
suggesting  that  particles  moving  with  the  group  velocity  can  lead 
to  break  up  and  dissipation  of  nonlinear  wavepackets.   The  matrix 
element  for  nonlinear  Landau  damping  extracted  from  this  theory  is 
shown  to  agree  in  the  appropriate  limit,  with  that  of  previous 
theories . 

Mr 

Work  supported  by  the  Center  for  Theoretical  Physics,  University 
of  Maryland,  and  by  the  Air  Force  Office  of  Scientific  Research. 
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Shock  Wave  Structure  in  a  Two  Component  Plasma 
Ferdinand  V".  Coronltl,  Burton  D.  Fried  and  Roscoe  B.  White 
University  of  California,  Los  Angeles 

Nonlinear  fluid  theory  equations  have  been  used  In  many 
instances  to  discuss  the  structure  of  the  leading  edge  of  a  shock 

wave,  or  In  the  absence  of  dissipation,  a  solitary  pulse.    Recently, 

2 
the  use  of  a  light  Ion  tracer  as  a  diagnostic  has  been  reported  and 

offers  the  possibility  of  a  detailed  examination  of  the  mechanisms 

Involved  in  collislonless  shock  structure.   In  this  paper  we 

investigate  the  effect  of  the  addition  of  a  light  ion  contaminant 

on  the  previous  treatments  of  shock  structure  referred  to  in  (1). 

The  nonlinear  coupled  fluid  and  kinetic  equations  are  solved, 

the  light  ions  treated  kinetically  and  the  heavy  ions  as  a  fluid,  and 

the  resulting  shock  structure  discussed  as  a  function  of  T  /T.,  mass 

ratio,  and  contamination.   The  two-fluid  limit  of  cold  light  ions 

is  also  obtained.   The  analysis  leads  to  a  new  critical  Mach  number, 

M  as  a  function  of  the  percentage  light  ion  contamination,  the 

ratio  T  /T.,  and  the  mass  ratio  of  the  two  ion  species.   It  describes 

the  transition  between  light  ion  reflection  by  the  shock  front  and 

the  penetration  of  the  front  by  the  light  ions.   Significant  changes 

in  shock  amplitude  and  leading  edge  profile  occur  at  this  critical 

mach  n^umber,  particularly  at  large  values  of  T  /T.  .   Shocks  of 

greatly  differing  amplitudes  (as  much  as  a  factor  of  10  for  an  A-He 

plasma)  may  be  found  for  the  same  values  of  light  ion  contamination 

and  mach  number  by  varying  T  /T. . 


Moiseev  and  Sagdeev,  Plasma  Phys .  _5,  ^3  (63)' 
Sagdeev,  Rev.  Plasma  Phys.  IV,  23  (66). 

2 

K.  R.  MacKenzie  and  R.  J.  Taylor,  Amer.  Phys.  Soc.  Meeting, 

November,  1970,  paper  2B4. 
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Evolution  of  Unstable  Bernstein  Modes  in  Astron  Plasma 
Charles  D.  Striffler  and  Terry  Kammash 
University  of  Michigan,  Ann  Arbor,  Michigan 
and  Marvin  Rensink 
LRL,  Livermore,  California 

Using  the  linearized  Vlasov  equations  for  the  Astron  system  we 
investigate  the  stability  of  the  Bernstein  modes  of  the  background 
plasma  as  a  result  of  interaction  with  the  relativistic  E- Layer,   We 
find  that  this  interaction  occurs  at  those  cyclotron  harmonics  of 
the  beam  which  lie  in  a  range  above  each  harmonic  of  the  background 
electron  cyclotron  frequency.   For  the  Bernstein  modes  between  the 
first  and  second  electron  harmonics  the  interaction  occurs  in  the 
frequency  range  between  the  electron  cyclotron  frequency  and  the 
upper  hybrid  frequency.   Considering  the  system  to  evolve  from  zero 
plasma  density  the  critical  density  necessary  for  the  onset  of  an 
instability  at  a  fixed  beam  harmonic  is  calculated.   The  growth  rate 
is  also  calculated  as  a  function  of  the  density  and  it  is  shown  that 
at  each  beam  harmonic  the  growth  rate  decreases  rapidly  with 

increasing  background  density.   These  results  agree  favorably  with 

2 
recent  experimental  observations  in  Astron  . 


Work  supported  by  the  U.  S.  Atomic  Energy  Commission 
■"■I.  B.  Bernstein,  Phys.  Rev.  109,  10  (1958). 
^T.  J.  Fessenden  and  B.  W.  Stallard,  UCRL-50002-70  (1970) 
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A  Schrodinger  Description  of  Non-Adiabatic  Particle  Loss 
and  of  Magnetic  Surfaces 
R.  K.  Varraa  and  C.  ¥.  Horton,  Jr. 
University  of  Texas  at  Austin 

A  model  is  obtained  which  describes  the  non-adlabatic  loss 
of  particles  from  magnetic  mirror  traps  through  Schrodinger- like 
equations 

#^=-^  ilif^^^m    ^(n),  n  =  1,2,3,... 

ox 

with  the  probability  density  G(x, t)  being  given  by 


G(x,t)  =XZ  ^  (n)^(n)  ; 
n 


here  the  role  of  n  is  played  by  the  first  action  invariant  il. 
The  process  of  non-adlabatic  loss  is  thus  found  to  be  analogous 
to  the  tunnelling  of  potential  barriers  in  quantum  mechanics,  the 
potential  (il  fi)  here  being  the  potential  for  the  adiabatlc  motion. 
The  predictions  of  this  model  for  the  life  times  compare  very  well 
with  experimental  results.   An  application  of  these  methods  for 
the  study  of  the  structure  of  magnetic  surfaces  is  indicated. 
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Decay  of  Pair  Correlations  in  a  Classical  Fluid 

J.  K.  Percus 
Courant  Institute  of  Mathematical  Sciences 
G.  J.  Yevick 
Stevens  Institute  of  Technology 

The  spatial  Fourier  components  of  the  pair  correlation 
function  of  a  classical  fluid  oscillate  in  time  over  a  considerable 
range  of  wave  vector  and  time,  lending  credence  to  an  analysis  in 
terms  of  appropriate  normal  modes.   The  decay  of  these  modes  can 
be  described  heuristically  by  a  phase  mixing  of  the  corresponding 
oscillators,  and  results  in  an  anomalous  behavior  for  plasmas, 
from  the  point  of  view  of  simple  Vlasov  theory.   To  approach  the 
problem  systematically,  it  is  convenient  to  focus  attention  upon  the 
particle  density  and  upon  a  suitably  defined  microscopic  velocity 
potential.   A  straightforward  singular  perturbation  theory  about 
the  independent  modes  can  then  be  carried  out  to  ellicit  the  long 
time  behavior.   However,  the  equal  time  values   of  the  velocity 
potential  correlations  are  required  as  input  data.   It  is  pointed 
out  that  the  assumption  of  universal  Gaussian  correlation  both 
reproduces  the  leading  order  perturbation  theory  results  and 
allows  direct  computation  of  the  relevant  equal  time  correlations. 
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Neoclassical  Electrical  Conductivity  of  Plasmas  with  Trapped  Particles 

Umur  Daybelge 
Physics  Department,  Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Plasmas  in  strong,  spatially  inhomogeneous  magnetic  fields 
have  a  lower  electrical  conductivity  than  plasmas  in  homogeneous 
fields.   For  such  plasmas  a  solution  method  for  the  kinetic 
equation  with  Fokker- Planck  operators  is  given.   The  non-local 
electrical  conductivity  is  calculated  and  compared  with  previous 
results . 
■X-  ... 

Visiting  Scientist 
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Self  Consistent  Field  Theory  of  Relativistic  Electron  Rings 

George  Schmidt 
Stevens  Institute  of  Technology 

Self  consistent  compression  of  a  thin  relativistic  electron 
ring  by  an  azimuthally  symmetric  confining  field  is  considered. 
The  relativistic  particle  Hamiltonian  is  expressed  in  terms  of  an 
effective  potential  in  two  dimensions,  a  function  of  the  constant 
canonical  angular  momentum  and  the  vector  and  scalar  potentials 
which  are  determined  self  consistently. 

For  a  neutralized  beam  the  self  consistent  ring  radius  is 
found  to  be  smaller  than  the  low  P  value  (contrary  to  naive 
intuition),  while  for  a  pure  electron  ring  the  radius  expands.   The 
particle  energy  always  decreases  due  to  collective  effects.   The 
variation  of  the  ring  cross  section  during  compression  is  estimated 
via  an  adabatic  invariant  associated  with  the  oscillating  particle 
in  the  effective  potential  well.   Curves  describing  the  variation 
of  these  quantities  in  terms  of  beam  and  field  parameters  have 
been  obtained  numerically.   Implications  to  possible  fusion  devices 
are  discussed. 
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Magnetically  Insulated  Ultra-High  Voltage  Transformer  and 
Controlled  Fusion  Research 
F.  Winterberg 
University  of  Nevada  System,  Las  Vegas,  Nevada 

If  the  magnetic  field  insulation  principle   is  applied  to  a 

2 
transformer  very  high  voltages  seem  to  be  attainable  .   With  the 

aid  of  such  a  transformer  continuous  intense  streams  of  relativistic 

electron  beams  could  be  extracted  from  the  high  voltage  terminal 

by  the  field  emission  process.   These  electron  streams  could 

be  efficiently  used  for  heating  a  plasma  confined  in  a  magnetic 

field,  especially  for  solving  the  energy  injection  problem  in 

open  ended  systems.   In  this  case  the  thermonuclear  fuel  injection 

would  have  to  be  simultaneously  accomplished,  which  could  be  done 

most  expediently  by  low-energy  large-mass  injection  with  a  catapult, 

propelled  by  a  superconducting  magnetic  travelling  wave  macron 

accelerator-^. 

The  electron  streams  produced  with  such  a  transformer  could 

be  also  used  to  charge  up  the  proposed  superconducting  levitated 

ring  capacitor  ,  making  possible  the  production  of  intense 

relativistic  electron  beams  with  energies  exceeding  anything  so 

far  achieved  and  thus  raising  the  prospect  of  thermonuclear 

1  4 
micro-bomb  ignition  '  . 

■"■F.  Winterberg,  Phys .  Rev.  ljj}_,    212  (1968). 

^F.  Winterberg,  Rev.  Scientific  Instr.  4l,  1756  (197O). 

^F.  Winterberg,  Plasma  Phys.  8,  5^1  (1966) . 

F.  Winterberg,  Proceedings  of  the  International  School  of  Physics 
"Enrico  Fermi"  on  High  Ene rgy  Density,  Varenna,  1969,    to  be 
published  by  the  Italian  Physical  Society. 
Supported  in  part  by  the  U.  S.  Air  Force  under  grant  No.  AFOSR-71-1980 
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Numerical  Simulation  of  Strong  Plasma  Shock  Waves 
Produced  In  An  Electromagnetic  Shock  Tube 
S.  H.  Schneider,  C.  K.  Chu 
Columbia  University 

A  set  of  two-fluid  Navier-Stokes  equations  with  classical 
physical  transport  coefficients  is  used  to  compute  the  evolution 
and  structure  of  collisional  plasma  shock  waves  in  an  electro- 
magnetic shock  tube.   Shock  speeds  up  to  ~  200  cm/iJ.  sec  and 
shocked  plasma  temperatures  of  the  order  of  kilovolts  are 
studied.   A  strong  transverse  bias  magnetic  field  is  employed, 
which  significantly  alters  the  size  and  shape  of  the  shock 
profiles,  when  compared  to  the  zero  bias  field  case.   The  wave 
structure  is  different  from  that  of  a  two-fluid  gas  dynamic 
shock  without  transverse  magnetic  field.   Especially  significant 
is  the  effect  of  the  small  ion  Larmor  radius  in  reducing  the 
gas  dynamic  shock  thickness  by  at  least  an  order  of  magnitude  in 
the  transverse  shock  case,  which  permits  collisional  shocks  to 
have  thicknesses  much  smaller  than  the  post-shock  mean  free  path. 
These  collisional  shock  waves  produce  very  hot  ions,  the  ion 
temperature  increasing  with  the  shock  speed  and  the  formation 
time  and  distances  also  increasing  substantially  with  shock 
speed.   Generally  the  results  compare  favorably  with  electro- 
magnetic shock  tube  experimental  data. 


Supported  by  the  United  States  Atomic  Energy  Commission  under 


Contract  AT( 30-1) 395^. 
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The  Role  of  Trapped  Particles  in  Plasma  Instabilities 

B.  Rosen,  G.  Schmidt 
Stevens  Institute  of  Technology 
W.  L.  Kruer 
Princeton  University 

We  report  on  a  number  of  Numerical  experiments  dealing 
with  the  role  that  trapped  particles  (TP)  play  in  the  saturation, 
damping  and  generation  of  plasma  instabilities.   It  is  found 
that  returning  the  TP  to  the  main  body  of  the  distribution 
eliminates  amplitude  oscillations  and  side  band  growth  for 
both  Landau  damping  and  beam  plasma  instability,   A  redistribution 
of  TP  over  the  trapping  region  of  phase  space  leads  to  significant 
changes  in  damping  rates  and  saturation  levels,  and  can  also 
alter  the  phase  of  the  amplitude  oscillation. 


Work  supported  in  part  by  the  AEG  and  NSF. 
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Low  Frequency  Heating 
William  B.  Thompson 
University  of  California  at  San  Diego 
La  Jolla,  California 

Several  problems  relevant  to  the  low  frequency  heating 
of  devices,  such  as  doublet,  having  high  magnetic  shear  will 
be  considered. 

* 
Supported  by  the  United  States  Atomic  Energy  Commission,  Contract 

No.  AT(04-3)-34  PA  85-13. 
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Theory  and  Simulation  of  the  Beam  Cyclotron  Instability 
Jl.    Lampe,  J.  B.  McBride,  W.  M.  Manheimer,  J.  H.  Orens,  R.  N.  Sudan  ,  R.  Shanny 

Naval  Research  Laboratory 

The  relative  drift  of  electrons  and  ions  across  a  magnetic 

1-4 
field  gives  rise  to  electron  cyclotron  instabilities.     We  present 

the  results  of  computer  experiments  for  such  modes  and  an  analysis 
which  explains  the  salient  features  of  the  observations.   Exponential 
quasilinear  heating  follows  the  onset  of  instability.   An  electro- 
static energy  in  the  wave  fields  builds  up,  the  electron  cyclotron 
resonances  are  broadened  by  turbulent  scattering.    At  a  level  of 

turbulence  y~"  I E,  1 2/4  ttN  T  (0)  ~  {^^1/^   ^^)[T  (t)/T  (0)]^/^(n  /kV,  .^.. 
*-7— '  k'     o  e^  ^    ^  ce  pe  ''-  e^  '   e^  ^     ^  ce   drift), 

where  T  (t)  is  the  electron  temperature  at  time  t,  the  effective 
collision  frequency  v  becomes  comparable  to  ^  ,  and  a  transition  is 
made  from  the  beam  cyclotron  instability  to  the  ion  acoustic  insta- 
bility.  The  system  then  continues  to  evolve  exponentially,  but  now 
at  the  slower  acoustic  rate,  until  the  fluctuating  fields  become 
strong  enough  to  trap  the  ions  and  stabilize  the  instability.   The 

final  temperatures  of  the  electrons  and  ions  are  given,  e.g., 

~       2   2 
T   <  (5/l6)  MV,  .^,,  and  the  important  consequences  for  perpendicular 

magnetosonic  shocks  of  the  above  behavior  are  discussed. 


On  leave  of  absence  from  Cornell  University. 

-•-H.  V.  Wong,  Phys.  Fluids  33,  757  (1970). 

2 
D.  W.  Forslund,  R.  L.  Morse  and  C.  W.  Niels on,  Phys.  Rev. 

Lett.  25,  1266  (1970) . 

-^S.  P.  Gary  and  J.  J.  Sanderson,  J.  Plasma  Phys.  4,  739  (1970)  . 

4 
M.  Lampe,  J.  B.  McBride,  J.  H.  Orens,  and  R.  N.  Sudan,  Phys. 

Lett.  A  (to  be  published). 

15 

-^We  use  the  resonance  broadening  equations  of  C.  T.  Dam  and  T.  H. 

Dupree,  Phys.  Fluids  13,  2064  (1970). 
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Electron  Vortices  and  Anomalous  DC  Resistivity 
W.  L.  Kruer  and  J.  M.  Dawson 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  N.  J.  085^0 

B.  Rosen 
Stevens  Institute  of  Technology,  Hoboken,  N.  J.  O703O 

In  computer  simulation  we  have  found  that  a  bifurcated 
electron  velocity  distribution  (two-symmetrical  electron  beams) 
results  in  a  substantial  transfer  of  energy  to  the  ions.   The 
energy  transfer  is  a  very  weak  function  of  the  electron  ion  mass 
ratio.   In  simulations  employing  realistic  electron-ion  mass 
ratios  we  have  observed  a  transfer  to  the  ions  of  up  to  13%   of 
the  total  electron  beam  energy.   A  simple  theory  explains  these 
results  in  terms  of  electron  vortices  in  phase  space  created  by 
the  strong  electron-electron  instability.   These  vortices  give  rise 
to  a  number  of  other  interesting  phenomena,  including  an  anomalous 
dc  resistivity.   When  we  apply  a  small  dc  field  to  this  system,  we 
observe  an  anomalous  resistance  for  the  lifetime  of  the  vortices. 
The  current  for  a  time  hangs  up  at  a  value  much  less  than  that 
determined  by  the  electron  thermal  velocity  (this  value  was  found 
in  previous  simulations  of  plasmas  driven  by  dc  fields).   We 
present  a  simple  picture  of  how  vortices  can  effectively  cause 
this  behavior. 


J.  P.  Boris,  J.  M.  Dawson,  J.  H.  Orens,  and  K.  V.  Roberts,  Phys. 
Rev,  Letters  25,  706  (1970). 

This  work  vjas  supported  jointly  by  U.  S.  Atomic  Energy  Commission 
Contract  AT(30-1)-1238  and  Office  of  Naval  Research  Laboratory 
Contract  N  00014-67-A-0151-0021. 
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Relatlvistic  Particle  Motion  in  Super- Intense 
Laser  Beams 
P,  K.  Kaw  and  R.  M,  Kulsrud 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  N.  J,  085^0 

Recently,  it  has  been  proposed  that  cosmic  rays  can  be 
accelerated  to  high  energies  in  the  intense  low-frequency  electro- 
magnetic waves  radiated  by  pulsars.   The  condition  for  relativistic 
acceleration  is  a)„  >  03  where  co  =  eB/mc  is  the  cyclotron  frequency 

rl  ri 

in  magnetic  field  B  of  the  wave  and  co  is  the  wave  frequency.   The 
same  conditions  can  be  achieved  in  the  laboratory  by  focusing  laser 
beams  with  powers  of  order  10   watts,  or  more,  down  to  a  few  wave 
lengths . 

We  have  examined  charged  particle  motions  in  focussed  laser 
pulses .   The  particle  motion  can  be  separated  into  a  guiding  center 
motion  and  a  figure  eight  gyration  analogous  in  many  ways  to 
particle  motion  in  a  magnetic  field.   Using  three  results,  we 
discuss  conditions  under  which  effective  acceleration  of  MeV 
electrons  can  be  achieved  in  the  laboratory.   Such  an  acceleration 
may  have  possible  application  to  Astron. 


This  work  was  supported  by  U.  S.  Air  Force  Office  of  Scientific 


Research  Contract  F44620-70-C-0033. 
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Toroidal  Effects  on  Kink  Modes  In  Tokamaks 
E.  A.  Frieman,  J.  M.  Greene,  and  K.  E.  Weimer 
Plasma  Physics  Laboratory,  Princeton  University 
Princeton,  New  Jersey  0854o 

Considerations  of  future  proposed  tokamaks  point  to  the 
desirability  of  achieving  smaller  aspect  ratios  R/a  and  lower 
values  of  the  safety  factor  q.   Previous  investigations  of 
tokamak  kink  stability  have  been  carried  out  in  cylindrical 
geometry,  thus  neglecting  the  magnetic  well  effects  of  toroidal 
geometry  and  the  possible  coupling  of  modes  of  differing  m 

numbers.   The  present  calculation  is  based  on  the  tokamak  equilib- 

2 
rium  expansion  of  Greene,  Johnson,  and  Weimer,   which  has  been 

carried  out  to  third  order  a/R.   To  go  beyond  the  cylindrical 

stability  theory  requires  an  expansion  to  fourth  order  in  a/R  as 

in  the  results  of  Ware  and  Haas-^  who  examined  tokamak  flute 

stability.   The  present  calculations  indicate  the  relation  between 

Pg,  q,  a/R,  and  various  higher  order  equilibrium  properties  for 

pressure  profiles  which  differ  in  higher  order  from  parabolic.   The 

results  are  obtained  by  assuming  the  second  order  neutrality 

requirements  are  satisfied  and  then  finding  corrections  to  these 

criteria. 


■'■v.  D.  Shafranov,  Zh.  Tekh.  Fiz.  4o,  24l  (197O)  [Sov.  Phys . 

Tech.  Phys.  15,  175(1970)]. 

2 
J.  M.  Greene,  J.  L.  Johnson,  and  K.  E.  Weimer,  Phys.  Fluids 

l4,  671  (1971). 

^A.  A.  Ware  and  F.  A.  Haas,  Phys.  Fluids  9,    956(1966). 

This  work  was  supported  by  the  U.  S.  Atomic  Energy  Commission, 

Contract  No.  AT( 30-l)-1238. 
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Computer  Simulation  of  Drift-Cyclotron  Instability 
H.  Okuda,  B.  Rosen,  and  J.  M,  Dawson 
Plasma  Physics  Laboratory,  Princeton  University,  Princeton,  N.J.  085^0 

A  study  of  the  drift-cyclotron  instability  and  the  associated 
plasma  diffusion  by  computer  simulation  has  been  started.   The  model 
used  is  electrostatic,  two  and  a  half  dimensional,  two  species  with 
uniform  external  magnetic  field.   The  code  is  optimized  and  with 
16,500  particles  of  each  species  on  a  32  x  32  grid  takes  0.4  8  seconds 
to  complete  one  time  step.   The  plasma  has  an  initial  density  gradient 
in  the  x  direction  with  the  guiding  center  distribution  ~  exp(-Kx)  and 
is  uniform  in  y  direction.   The  velocity  distribution  of  the  ions 
and  electrons  are  Maxwellian  with  the  same  temperature.   This  differs 
from  the  recent  results  by  Birdsall  et  al,   who  examined  the  drift 
cyclotron  instability  for  monoenergetic  ions  and  cold  electrons. 
When  k||  =  0,  an  instability  arises  due  to  the  coupling  between 

negative  energy  ion  cyclotron  harmonics  and  the  electron  drift  wave 

2 
with  a  maximum  growth  rate  of  y/(a   .    «  0.5  for  (co  ./co  .)   =  25,  of 
&  '  ci  ^  pi  ci'     -^' 

Kp  2;  1,  and  m./m  =4.   All  unstable  modes  grow  with  growth  rates 
close  to  the  predicted  value.   The  most  unstable  mode  reaches  a 
maximum  amplitude  of  e<!>/kT  x  4o  at  co  .  t  x  30,  followed  by  oscillations 
of  its  amplitude.   With  the  development  of  large  waves,  electrons  and 
then  ions  diffuse  and  spread  out  in  the  x  direction.   However,  the 
large-amplitude  waves  in  the  y  direction  cause  a  large  density 
modulation  in  that  direction  producing  an  electric  field  which  in 
turn  causes  electrons  to  drift  in  x  direction  and  to  some  extent 
restores  the  original  density  profile. 


-'"C.  K.  Birdsall  and  D.  Fuss,  Am.  Phys .  Soc.  4C12,  1437(1970). 

This  work  was  supported  by  U.  S.  Atomic  Energy  Commission  Contract 

AT( 30-I) -1238;  use  was  made  of  computer  facilities  supported  in  part 

by  National  Science  Foundation  Grant  NSF-GP  579. 
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Analysis  of  Nonlinear  Plasma  Waves 
James  T.  Yen 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

and 
Grumman  Research  Laboratories,  Bethpage,  N.Y. 

The  exact  solution  to  the  Vlasov  equation  can  be  expressed 
in  terms  of  an  exponential  operator,  operating  on  the  initial 
velocity  distribution.   If  we  take  small  time-steps,  this 
exponential  operator  can  be  approximated  by  a  non-linear  poly- 
nomial operator.   The  latter  is  used  to  obtain  analytically  at 
subsequent  times,  the  velocity  distribution,  the  electric  field, 
the  non-linear  Landau  damping  or  amplification,  and  the  diffusion 
of  the  plasma.   The  plasma  will  occupy  a  finite  one-dimensional 
region  with  assigned  boundary  conditions;  the  motion  of  the 
ions  will  be  neglected. 

Starting  with  a  single  wave  number  and  initial  velocity 
distribution  containing  a  cubic  power  small  disturbance,  the 
analysis  readily  produces  all  wave  numbers.   By  expressing 
non-linear  exponential  functions  in  terms  of  infinite  series  of 
Bessel  and  circular  functions,  we  can  analyze  exactly  the  non- 
linear mode-mode  coupling  among  all  the  wave  numbers.   The 
boundary  conditions  will  be  brought  in  exactly.   The  change  in 
shape  and  width  of  the  disturbance  will  be  demonstrated. 
Diffusion  and  nonlinear  damping  or  amplification  is  studied. 
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Energy  Principle  with  Given  Distributed  Inductance 

Glenn  Bateman 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

Holding  the  shape  of  a  set  of  nested  flux  surfaces  fixed, 
we  study  the  variation  of  the  ideal  MHD  potential  energy  with 
respect  to  rearrangements  of  the  flux  and  pressure  profiles. 
The  constraint  on  the  shape  of  the  flux  surfaces  is  incorporated 

r    12 

by  writing  /  drpB  as  an  integral  over  a  bilinear  form  of  flux 
and  the  inverse  of  the  differential  inductance  matrix.   The 
inductance  matrix  is  held  fixed  while  the  flux  profile  is  varied, 
A  number  of  thermodynamic  models  are  used  in  the  variation  of 
the  pressure  profile;  each  model  gives  a  necessary  stability 
condition.   For  the  models  considered,  and  for  a  diagonal  in- 
ductance matrix,  the  state  of  minimum  energy  is  a  constant- 
pressure  sharp  boundary  plasma  with  only  longitudinal  current. 
An  azimuthal  current  increases  stability  with  respect  to  per- 
turbations in  the  flux  profile;  smoothing  the  pressure  profile 
can  decrease  stability. 
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Hydromagnetic  Stability  of  Bifurcated  Equilibria 

Tyan  Yeh 
Courant  Institute  of  Mathematical  Sciences,  New  York  University 

Helically  sjaranetric  hydromagnetic  equilibria  are  considered 
in  a  cylindrical  domain.   These  non-unique  solutions  may  be 
regarded  as  bifurcated  from  a  straight  pinch.   Such  a  bifurcated 
equilibrium  is  generally  unstable  unless  the  bifurcation  point 
occurs  at  the  margin  of  stability,  in  other  words,  at  a  higher 
order  confluence  of  bifurcating  equilibria.   In  the  latter  case, 
we  show  that  the  bifurcating  equilibrium  is  stable  when  the 
sharp  plasma  vacuum  interface  has  a  long  helical  wavelength. 

First  we  expand  the  equilibrium  solution  in  a  small 
parameter,  the  amplitude  of  the  corrugation  of  the  plasma-vacuum 
interface.   Then  we  identify  this  class  of  helically  symmetric 
equilibria  as  bifurcated  equilibria  and  discuss  a  necessary 
condition  for  the  bifurcated  equilibrium  to  be  stable.   Finally 
using  the  energy  principle  for  hydromagnetic  stability  we  study 
the  stability  of  the  bifurcated  equilibria. 

It  is  an  open  question  whether  the  new  class  of  bifurcated 
equilibria  is  stable  when  the  amplitude  of  the  corrugation  is 
large  and  when  the  helical  wavelength  is  not  large. 
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This  report  was  prepared  as  an  account  of 
Government  sponsored  work.   Neither  the 
United  States,  nor  the  Commission,  nor  any 
person  acting  on  behalf  of  the  Commission: 

A.  Makes  any  warranty  or  representation, 
express  or  implied,  with  respect  to  the 
accuracy,  completeness,  or  usefulness  of 
the  information  contained  In  this  report, 
or  that  the  use  of  any  information, 
apparatus,  method,  or  process  disclosed 
In  this  report  may  not  infringe  privately 
owned  rights;  or 

B.  Assumes  any  liabilities  with  respect  to 
the  use  of,  or  for  damages  resulting  from 
the  use  of  any  information,  apparatus, 
method,  or  process  disclosed  in  this 
report. 

As  used  in  the  above,  "person  acting  on  behalf 
of  the  Commission"  Includes  any  employee  or 
contractor  of  the  Commission,  or  employee  of 
such  contractor,  to  the  extent  that  such  em- 
ployee or  contractor  of  the  Commission,  or 
employee  of  such  contractor  prepares,  dis- 
seminates, or  provides  access  to,  any  infor- 
mation pursuant  to  his  employment  or  contract 
with  the  Commission,  or  his  employment  with 
such  contractor. 
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